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ABSTRACT 


This  thesis  investigates  the  potential  use  of  wurtzite 
Indium  Gallium  Nitride  as  photovoltaic  material.  Silvaco 
Atlas  was  used  to  simulate  a  quad- j unction  solar  cell.  Each 
of  the  junctions  was  made  up  of  Indium  Gallium  Nitride.  The 
band  gap  of  each  junction  was  dependent  on  the  composition 
percentage  of  Indium  Nitride  and  Gallium  Nitride  within 
Indium  Gallium  Nitride.  The  findings  of  this  research  show 
that  Indium  Gallium  Nitride  is  a  promising  semiconductor  for 
solar  cell  use. 
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EXECUTIVE  SUMMARY 


One  of  the  primary  goals  of  solar  cell  design  is  to 
improve  efficiency.  Indium  Gallium  Nitride  is  a  material 
that  has  undergone  extensive  research  since  2002  as  a 
potential  photovoltaic  material.  By  varying  the  composition 
of  Indium  Nitride  and  Gallium  Nitride  within  Indium  Gallium 
Nitride,  the  band  gap  of  this  semiconductor  material  can  be 
changed.  The  band  gap  range  of  Indium  Gallium  Nitride 
matches  closely  the  visible  solar  spectrum  frequencies. 
Hence,  a  high-efficiency  solar  cell  can  be  potentially 
developed  by  having  several  Indium  Gallium  Nitride 
junctions . 

The  use  of  Silvaco  Atlas  as  a  simulation  tool  can  aid 
in  identifying  whether  Indium  Gallium  Nitride  can  be  used 
for  photovoltaics .  Previous  research  at  the  Naval 
Postgraduate  School  has  demonstrated  the  viability  of 
Silvaco  Atlas  as  a  solar  cell  modeling  mechanism. 

Finding  high-efficiency  solar  cell  models  is  of  great 
interest  in  space  applications.  By  increasing  the  efficiency 
of  photovoltaics,  the  number  of  solar  panels  is  decreased. 
Therefore,  the  overall  weight  that  needs  to  be  launched  into 
space  is  reduced.  A  cost-reduction  can  be  achieved  over  the 
life  of  the  space  power  system. 

Research  at  Lawrence  Berkeley  National  Laboratory 
(LBNL)  has  been  ongoing  since  2002  in  trying  to  develop 
high-efficiency  solar  cells.  One  of  the  materials  LBNL  has 
been  investigating  is  Indium  Gallium  Nitride.  In  addition  to 
the  band  gap  range.  Indium  Gallium  Nitride  also  has  other 
characteristics  that  are  beneficial  for  space  applications. 
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The  radiation  level  that  Indium  Gallium  Nitride  is  able  to 
withstand  is  approximately  two  orders  of  magnitude  greater 
than  current  multi j unction  solar  cell  materials.  Therefore, 
it  is  of  great  interest  to  continue  to  research  Indium 
Gallium  Nitride  and  its  possibilities  as  photovoltaic 
material . 

The  goal  of  this  research  was  to  investigate  the 
efficiency  of  Indium  Gallium  Nitride  solar  cells  using  the 
Silvaco  Atlas  TCAD  simulation  software.  No  other  TCAD 
simulations  have  been  performed  on  InGaN  solar  cells.  The 
methodology  of  this  research  progressed  from  single  to 
multi- j unction  cells.  The  simulation  results  predicted 
efficiencies  as  high  as  41%  for  a  four- j unction  solar  cell. 
Further  refinements  of  the  simulation  model  are  still 
possible.  Actual  production  of  a  single- j unction  solar  cell 
is  the  next  step  required  to  eventually  manufacture  high- 
efficiency,  multi- j unction  InGaN  solar  cells. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

One  of  the  most  important  problems  to  solve  in  space 
solar  cell  design  is  efficiency.  As  efficiency  increases, 
the  required  number  of  cells  decreases  to  fulfill  electrical 
power  requirements. 

Solar  cell  efficiencies  have  improved  over  time  by 
increasing  the  number  of  junctions.  Each  junction  is  capable 
of  extracting  energy  from  a  portion  of  the  solar  spectrum. 

A  new  path  to  improving  efficiency  is  to  use  new 
photovoltaic  materials  such  as  wurtzite  Indium  Gallium 
Nitride.  Indium  Gallium  Nitride  possesses  a  band  gap  range 
that  can  extract  energy  from  a  large  portion  of  the  visible 
solar  spectrum. 

B.  OBJECTIVE 

The  objective  of  this  thesis  is  to  simulate  Indium 
Gallium  Nitride  solar  cells  to  predict  their  efficiencies. 
After  obtaining  the  efficiencies,  a  comparison  is  made  with 
current  solar  cell  efficiencies.  Silvaco  Atlas  TCAD 
simulation  software  is  used  to  investigate  these 
efficiencies . 

C .  RELATED  WORK 

Previous  Silvaco  Atlas  solar  cell  simulations  have  been 
performed  by  Naval  Postgraduate  School  researchers. 
Michalopoulos  [1]  investigated  the  feasibility  of  designing 

solar  cells  using  Silvaco  Atlas.  To  demonstrate  the  use  of 
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Silvaco  Atlas,  Michalopoulos  simulated  single- j unction  solar 
cells  with  Gallium  Arsenide,  dual- j unction  solar  cells  with 
Indium  Gallium  Phosphide  and  Gallium  Arsenide,  and  triple¬ 
junction  cells  with  Indium  Gallium  Phosphide,  Gallium 
Arsenide,  and  Germanium.  The  highest  efficiency  obtained 
with  the  triple- j unction  was  29.5%.  This  result  matched  the 
29.3%  efficiency  obtained  with  actual  triple- j unction  solar 
cells  in  production.  Bates  [2]  provided  excellent  background 
on  the  use  of  Silvaco  Atlas.  Green  [3]  and  Canfield  [4] 
continued  to  work  on  Silvaco  Atlas  solar  cell  design.  Matlab 
functions  and  Silvaco  Atlas  code  were  modified  from  [  1 ]  —  [ 4 ] 
to  perform  simulations  in  this  research. 

A  previous  Indium  Gallium  Nitride  simulation  has  been 
performed  [5]  using  fundamental  semiconductor  physics 
formulas.  The  results  of  this  research  are  compared  with 
those  of  [ 5 ] . 

D.  ORGANIZATION 

1.  Purpose  of  Solar  Cells 

The  purpose  of  solar  cells  is  to  convert  light  energy 
into  electrical  energy.  Light  is  made  up  of  photons.  Photons 
carry  energy  that  is  dependent  on  the  color,  or  wavelength 
of  light.  Electrical  energy  is  generated  when  photons  excite 
electrons  from  the  valence  band  into  the  conduction  band  in 
semiconductor  materials.  Chapter  II  covers  fundamental 
concepts  of  solar  cells  and  semiconductor  materials. 

2.  Simulation  Software 

Silvaco  Atlas  was  used  as  the  simulation  software  in 

this  thesis.  Silvaco  offers  a  suite  of  software  programs 
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that  can  act  as  a  Virtual  Wafer  Fabrication  (VWF)  tool.  The 
Naval  Postgraduate  School  has  had  multiple  researchers 
publish  theses  on  the  use  of  Silvaco  Atlas  for  the  purpose 
of  solar  cell  modeling.  Chapter  III  covers  the  basics  of 
Silvaco  Atlas. 

3 .  Indium  Gallium  Nitride 

Wurtzite  Indium  Gallium  Nitride  (InGaN)  is  a 
semiconductor  that  has  the  potential  to  produce  high- 
efficiency  solar  cells.  Dr.  Wladek  Walukiewicz,  of  Lawrence 
Berkeley  National  Laboratory  (Solar  Energy  Materials 
Research  Group) ,  has  been  conducting  extensive  InGaN 
research  for  the  purposes  of  producing  photovoltaic 
material.  Chapter  IV  covers  the  characteristics  of  InGaN, 
its  potential,  its  material  science  research  status,  and  its 
current  limitations. 

4 .  Simulation 

Based  on  the  optical  characteristics  of  InGaN, 
simulations  were  performed  using  Silvaco  Atlas.  Chapter  V 
covers  the  findings  of  the  simulations.  Chapter  VI  covers 
the  conclusions  and  recommendations.  Finally,  the  appendices 
include  the  code  used  in  Silvaco  Atlas  as  well  as  Matlab. 
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II.  SOLAR  CELL  AND  SEMICONDUCTOR  BASICS 


This  thesis  examines  a  novel  process  to  simulate  Indium 
Gallium  Nitride  solar  cells.  Before  delving  into  the 
simulation,  this  chapter  covers  the  fundamental  principles 
of  solar  cells  and  semiconductors.  The  properties  of  the 
semiconductor  material  determine  the  characteristics  of  the 
photovoltaic  device. 

A.  SEMICONDUCTOR  FUNDAMENTALS 


1.  Classification  of  Materials 


Materials  can  be  categorized  according  to  their 
electrical  properties  as  conductors,  insulators  or 
semiconductors.  The  conductivity  cr  is  a  key  parameter  in 
identifying  the  type  of  material. 

Figure  1  presents  a  sample  of  materials  based  on 
conductivity.  The  semiconductors  fall  between  the  insulators 
and  the  conductors . 


Resistivity  p  (fl  -  cm) 
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Figure  1.  Materials  classified  by  conductivity  (From  [6]) . 
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Semiconductors  can  be  found  in  elemental  or  compound 
form.  Silicon  (Si)  and  Germanium  (Ge)  are  examples  of 
elemental  semiconductors.  Both  of  these  semiconductors 
belong  to  group  IV  of  the  periodic  table. 

Figure  2  shows  an  abbreviated  periodic  table.  In 
addition  to  the  group  IV  semiconductors,  compounds  can  be 
made  with  elements  from  groups  III  and  V,  respectively. 
Examples  of  III-V  semiconductors  include  Aluminum  Phosphide 
(A1P) ,  Gallium  Nitride  (GaN) ,  Indium  Phosphide  (InP), 
Gallium  Arsenide  (GaAs) ,  among  others.  It  is  also  possible 
to  make  semiconductor  compounds  from  groups  II-VI,  such  as 
Zinc  Oxide  (ZnO) ,  Cadmium  Telluride  (CdTe) ,  Mercury  Sulfide 
(HgS) ,  among  others.  The  subject  of  this  thesis  is  the 
ternary  alloy  wurtzite  Indium  Gallium  Nitride  (InGaN). 
Indium  and  Gallium  are  group  III  elements,  while  Nitrogen  is 
a  group  V  element. 


3A  | 

4A  | 

5A  [ 

6A  | 

5  B 

6  c 

7  N 

8  0 

10.811 

Boron 

12.0107 

Carton 
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15.9994 
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13  Al 

14  Si 

15  p 

16  § 

2bIH 

26.9815 

Aluminum 

28.0855 

Silicon 
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3°  zn 

31  Ga 

32  Ge 
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34  Se 
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48  Cd 

49  In 

50  Sn 

51  Sb 

52  Te 
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80  Hg 
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82  Rb 

83  Bi 

84 

200.59 
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207.2 
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208.9804 

Bismuth 

[209| 

Polonium 

Figure  2.  Partial  periodic  table  (After  [7]) . 
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2. 


Atomic  Structure 


Since  Silicon  is  the  most  commonly  used  semiconductor 
in  solar  cells  today,  a  brief  analysis  of  this  element  is 
presented . 

Silicon  has  14  protons  and  14  neutrons  in  its  nucleus. 
The  14  electrons  are  distributed  in  three  shells. 

Figure  3  shows  the  arrangement  of  electrons  in  a 
Silicon  atom.  There  are  two  electrons  in  the  first  shell, 
eight  electrons  in  the  second  shell  and  four  electrons  in 
the  outer  shell.  When  Silicon  atoms  are  together,  the  atoms 
from  the  outer  shells  form  covalent  bonds.  Hence,  a  Silicon 
atom  forms  bonds  with  four  other  Silicon  atoms. 


Figure  3.  Silicon  atomic  structure  (From  [8]) . 

As  shown  in  Figure  4,  the  covalent  bonds  formed  by  the 
silicon  atoms  are  represented  by  the  elliptical  dotted 
lines.  In  the  absence  of  an  electron,  the  covalent  bond 
ceases  to  exist.  A  hole  takes  the  place  of  the  electron. 
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Figure  4.  Silicon  atom  covalent  bonds  (From  [1,  p.  11]) . 

Energy  bands  are  a  fundamental  concept  of  semiconductor 
physics.  These  energy  bands  are  the  valence  band,  the 
conduction  band,  and  the  forbidden  gap  or  band  gap  [9,  p. 
27]  .  When  an  electron  is  in  the  valence  band,  the  covalent 
bond  exists.  In  order  for  the  electron  to  move  from  the 
valence  band  into  the  conduction  band,  energy  is  required  to 
excite  the  electron.  The  band  gap  energy  is  the  minimum 
energy  required  for  the  electron  to  make  the  move  from  the 
valence  band  into  the  conduction  band. 

Figure  5  shows  the  band  gap  diagrams  for  conductors, 
insulators,  and  semiconductors.  This  figure  also  reinforces 
the  differences  among  the  three  types  of  materials  according 
to  their  electric  properties.  In  the  case  of  conductors,  the 
band  gap  is  small  or  non-existent.  By  contrast,  insulators 
have  wide  band  gaps.  Therefore,  it  takes  much  more  energy 
for  the  insulator  to  have  electrons  in  the  conduction  band. 
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Figure  5.  Band  gap  diagrams  (From  [1,  p.  8]) . 

Semiconductors  have  band  gaps  that  are  dependent  on  the 
material.  Because  band  gap  is  also  dependent  on  temperature, 
it  should  be  noted  that  all  quoted  band  gaps  in  this  thesis 
are  at  room  temperature  (300  K.) 

3 .  Electrons  and  Holes 

A  semiconductor  at  absolute  zero  temperature  is  unable 
to  conduct  heat  or  electricity  [10,  p.  44]  .  All  of  the 
semiconductor's  electrons  are  bonded.  The  electrons  acquire 
kinetic  energy  as  the  temperature  is  increased.  Some  of  the 
electrons  are  freed  and  move  into  the  conduction  band.  These 
electrons  are  able  to  conduct  charge  or  energy.  The  areas 
left  by  these  electrons  in  the  valence  band  are  called 
holes.  As  the  temperature  of  the  semiconductor  increases, 
the  number  of  free  electrons  and  holes  increases  as  well. 
Therefore,  the  conductivity  of  the  semiconductor  is  directly 
proportional  to  temperature  increases. 

The  conductivity  of  the  semiconductor  can  be  altered  by 
exposing  it  to  light  or  by  doping  it.  Photoconductivity 
consists  of  exposing  the  semiconductor  with  photon  energy 
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larger  than  the  semiconductor's  band  gap.  Doping  consists  of 
adding  impurities  to  the  semiconductor  material. 


Figure  6  shows  n-type  and  p-type  doping.  The  atom  of 
the  semiconductor  is  represented  by  the  blue  circle  and  a 
+  4.  For  example,  a  silicon  atom  has  four  electrons  in  its 
outer  shell.  When  silicon  atoms  are  n-doped,  atoms  from 
group  V  of  the  periodic  table  are  added.  Each  added  group  V 
atom  provides  an  extra  electron  to  donate.  When  silicon 
atoms  are  p-doped,  the  atoms  from  group  III  of  the  periodic 
table  are  added.  Each  group  III  atom  has  one  less  electron 
than  silicon.  Therefore,  a  hole  is  added.  Doping  increases 
the  number  of  carriers  in  the  semiconductor  material. 
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Figure  6.  Doping:  n-type  and  p-type  (From  [1,  p.  13]). 


An  intrinsic  semiconductor  is  a  pure  semiconductor  with 
a  negligible  amount  of  impurity  atoms  [9,  p.34] .  By 
definition,  the  number  of  electrons  and  holes  in  a 
semiconductor  are  represented  by  n  and  p.  See  Table  1. 
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number  of  electrons 

n  = - 3 - 

cm 

number  of  holes 

P  = - I - 

_ cnr _ 

Table  1.  Definitions  of  n  and  p. 

In  the  case  of  an  intrinsic  semiconductor,  the 
following  case  occurs: 

n  =  p  =  m 

For  example  purposes,  ni  is  given  for  GaAs,  Si  and  Ge 
and  room  temperature  [9,  p.34]  .  Data  for  InN  and  GaN  are 

from  [11] .  See  Table  2. 


GaAs 

ni 

_  2xl06 

cm 3 

n? 

lxlO10 

cm 3 

nfe 

_  2xl013 

cnr’ 

njnN 

_  6.2xl022 

cm 3 

nGaN 

_  8.9 xlO22 

cm 3 

Table  2.  ni  for  five  semiconductors. 


When  solar  cells  utilize  semiconductor  materials  from 
Table  2,  the  current  level  is  highest  for  Ge  and  lowest  for 
GaAs.  The  current  levels  correspond  in  rank  to  the  n±  levels 
provided  in  Table  2.  No  data  exists  for  InN  or  GaN  current 
levels  as  photovoltaic  materials. 
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4. 


Direct  and  Indirect  Band  Gaps 


Since  the  band  gap  is  the  minimum  energy  required  to 
move  an  electron  from  the  valence  band  into  the  conduction 
band,  it  is  necessary  to  distinguish  between  direct  and 
indirect  band  gaps. 

Figure  7  shows  the  concept  of  direct  and  indirect  band 
gaps.  The  blue  portion  represents  the  valence  band.  The  tan 
portion  represents  the  conduction  band. 


Figure  7.  Direct  and  indirect  band  gaps  (After  [1,  p.  24]). 

When  the  valence  band  and  the  conduction  band  coincide 
in  wave  vector  k,  the  semiconductor  has  a  direct  band  gap. 
When  the  valence  band  and  the  conduction  band  have  different 
wave  vector  k,  the  semiconductor  has  an  indirect  band  gap. 
The  k  vector  represents  a  difference  in  momentum.  Photons 
have  negligible  momentum.  In  order  to  excite  an  electron 
from  the  valence  band  to  the  conduction  band  in  an  indirect 
band  gap  semiconductor,  in  addition  to  a  photon,  a  phonon  is 
required.  The  phonon  is  a  lattice  vibration.  The  phonon 


12 


transfers  its  momentum  to  the  electron  at  the  time  the 
photon  is  absorbed.  Therefore,  a  direct  band  gap 
semiconductor  is  generally  better  for  optoelectronics. 
Silicon  is  an  example  of  an  indirect  band  gap  semiconductor. 
Gallium  Arsenide  and  wurtzite  Indium  Gallium  Nitride  are 
examples  of  direct  band  gap  semiconductors. 

5 .  Fermi  Level 

The  Fermi  function  f (E)  specifies  how  many  of  the 
existing  states  at  energy  E  are  filled  with  an  electron  [9, 
p.42]  .  The  Fermi  function  is  a  probability  distribution 
function  defined  as  follows: 

1 

(. E-Ef ) 

1  +  e  kT 

Where  EF  is  the  Fermi  level,  k  is  Boltzmann's  constant, 
and  T  is  the  temperature  in  Kelvin. 

From  the  Fermi  function,  it  can  be  determined  that  when 
E=Ef,  then  f  (E)  =f  (EF)  =0 . 5  . 

Figure  8  shows  the  Fermi  level  for  an  intrinsic 

semiconductor.  Figure  9  shows  the  Fermi  level  for  n-type 
semiconductor.  Figure  10  shows  the  Fermi  level  for  p-type 
semiconductor.  From  these  Figures,  it  can  be  deduced  that 
the  n-type  material  has  a  larger  electron  carrier 

distribution  and  the  p-type  material  has  a  larger  hole 
carrier  distribution. 

- £c 

£F  . E\ 

- £v 
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Figure  8  . 


Fermi  level:  intrinsic  case  (After  [9,  p.  42]). 
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Figure  9.  Fermi  level:  n-type  case  (After  [9,  p.  42]). 

- Ec 

- E. 

- Ep 

- £v 

P-typc 

Figure  10.  Fermi  level:  p-type  case  (After  [9,  p.  42]). 

B.  SOLAR  CELL  FUNDAMENTALS 

After  covering  the  basics  of  semiconductors,  the 
logical  step  is  to  continue  with  solar  cell  fundamentals. 

1.  History  of  Solar  Cells 

Solar  cell  research  started  when  Edmund  Bequerel 
discovered  the  photovoltaic  effect  in  1839.  He  recorded  that 
an  electric  current  was  produced  when  light  was  applied  to  a 
silver  coated  platinum  electrode  immersed  in  electrolyte. 
The  next  significant  step  was  performed  by  William  Adams  and 
Richard  Day  in  1876.  They  discovered  that  a  photocurrent 
appeared  when  selenium  was  contacted  by  two  heated  platinum 
contacts.  However,  current  was  produced  spontaneously  by  the 
action  of  light.  Continuing  to  build  on  these  efforts, 
Charles  Fritts  developed  the  first  large  area  solar  cell  in 
1894.  He  pressed  a  layer  of  selenium  between  gold  and 
another  metal.  Progress  continued  as  the  theory  of  metal- 
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semiconductor  barrier  layers  was  established  by  Walter 
Schottky  and  other.  In  the  1950s,  silicon  was  used  for  solid 
state  electronics.  Silicon  p-n  junctions  were  used  to 
improve  on  the  performance  of  the  Schottky  barrier.  These 
silicon  junctions  had  better  rectifying  action  and 
photovoltaic  behavior.  In  1954,  Chapin,  Fuller  and  Pearson 
developed  the  first  silicon  solar  cell,  with  a  reported 
efficiency  of  6%.  However,  the  cost  per  Watt  associated  with 
these  solar  cells  made  them  prohibitively  expensive  for 
terrestrial  use.  However,  locations  where  power  generation 
was  not  feasible  (i.e.,  space)  were  suitable  for  solar 
cells.  Satellites  were  the  first  clear  application  for 
silicon  solar  cells.  Since  that  time,  solar  cells  have 
progressed  steadily  both  in  terms  of  efficiency  as  well  as 
materials  used  their  production  [10,  p.  2] . 

Table  3  shows  a  brief  list  of  events  in  the  history  of 
photovoltaics  from  1939  until  2002. 

2.  The  Photovoltaic  Effect 

The  photovoltaic  effect  is  the  process  by  which  a  solar 
cell  converts  the  energy  from  light  into  electrical  energy. 
Light  is  made  up  of  photons.  The  energy  of  these  photons 
depends  on  the  color  (wavelength)  of  light.  The  material 
that  makes  up  the  solar  cell  determines  the  photovoltaic 
properties  when  light  is  applied  [10,  p.  1] . 

When  light  is  absorbed  by  matter,  such  as  metal, 
photons  provide  the  energy  for  electrons  to  move  to  higher 
energy  states  within  the  material.  However,  the  excited 
electrons  return  to  their  original  energy  state.  In 
semiconductor  materials,  there  is  a  built-in  asymmetry  (band 
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gap)  .  This  allows  the  electrons  to  be  transferred  to  an 
external  circuit  before  they  can  return  to  their  original 
energy  state.  The  energy  of  the  excited  electrons  creates  a 
potential  difference.  This  electromotive  force  directs  the 
electrons  through  a  load  in  the  external  circuit  to  perform 
electrical  work. 


•  1839  Bevqucrel  <  I  K  i  discovered  photogalvanic  cITccl  in  liquid  electrolytes 

•  1873  Smith  (UKi  discovered  photoconductivity  of  solid  Se 

•  1877  Admits  and  Day  (UK)  discover  photogcncralion  of  eurrenl  in  Se  lubes:  the  lirsl  observation  of 
PV  effect  in  solids 

•  1883  Frills  (US)  makes  lirsl  large  area  solar  cell  using  Se  lilm 

•  1954  First  6r/J  efficient  solar  cells  reported:  Si  (Bell  Lab.  USA)  and  Cu:S/CdS  (Air  Force.  USA) 

•  1955  Hodman  Electronics  (USA)  offers  2'A  efficient  Si  PV  cells  at  S.I50OAV 

•  1958  NASA  Vanguard  satellite  with  Si  backup  solar  array 

•  1959  Hoffman  Electronics  (USA)  offers  105;  efficient  Si  PV  cells 

•  1963  Sharp  Corp  (JP)  produces  lirsl  commercial  Si  modules 

•  1966  NASA  Orbiting  Astronomical  Observatory  launched  with  I  k\V  array 

•  1970  First  GaAs  heterostructure  solar  cells  by  Alferov.  Andreev  el  at.  in  the  USSR 

•  1972  First  PV  conference  to  include  a  session  on  terrestrial  applications  (IEEE) 

•  1973  A  big  year  in  photovoltaics:  Worldwide  oil  crisis  spurs  many  nations  to  consider  renewable 
energy  including  photovoltaics:  Cherry  Hill  Conference  in  USA  (established  photovoltaics’  potential 
and  legitimacy  for  government  research  funding):  World's  lirsl  solar  powered  residence  (University  of 
Delaware.  USA)  built  with  Cu.S  (not  c-Si!)  solar  modules 

•  1974  Project  Sunshine  initialed  in  Japan  to  foster  growth  of  PV  industry  and  applications:  Tyco  (USA) 
grows  2.5  cm  wide  Si  ribbon  for  photovoltaics.  lirsl  alternative  to  Si  wafers 

•  1975  First  book  dedicated  to  PV  science  and  technology  by  Hovel  (USA) 

•  1980  First  thin-lilm  solar  cell  >  HI'/!  using  CU’S/CdS  (USA) 

•  1981  350  kW  Concentrator  array  installed  in  Saudi  Arabia 

•  1982  First  I  MW  utility  scale  PV  power  plant  (CA.  USA)  with  Arco  Si  modules  on  2-axis  trackers 

•  1984  6  MW  array  installed  in  Carrisa  Plains  CA.  USA  [ 35 1 

•  1985  A  big  year  lor  high-efficiency  Si  solar  cells:  Si  solar  cell  >  20'/!  under  standard  sunlight 
(UNSW.  Australia)  |36]  and  >25*4  under  200X  concentration  (Stanford  Univ.  USA)  |37| 

•  1986  First  commercial  thin-lilm  power  module,  the  a-Si  (i4(K)()  from  Arco  Solar  (USA) 

•  1987  Fourteen  solar  powered  cars  complete  the  3200  km  World  Solar  Challenge  race  (Australia)  with 
the  winner  averaging  70  kph 

•  1994  CialnP/CiaAs  2-terminal  concentrator  multifunction  >30';  (NREL.  USA)  |38| 

•  1995  "1000  roofs”  German  demonstration  project  to  install  photovoltaics  on  houses,  which  triggered 
the  present  favorable  PV  legislation  in  Germany.  Japan  and  other  countries 

•  1996  Photoelectrochemical  “dyc-sensiti/ed"  solid/liquid  cell  achieves  1 1'!  (EPFL.  Switzerland)  [39| 

•  1997  Worldwide  PV  production  reaches  100  MW  per  year 

•  1998  Cu(lnGa)Se.  thin-lilm  solar  cell  reaches  19',;  efficiency  (NREL.  US)  |40|  comparable  with 
mullicryslalline  Si.  First  concentrating  array  for  space  launched  on  Deep  Space  I  by  US  (5  kW  using 
high  efficiency  CalnP/GaAs/Gc  triple  junction  cells) 

•  1999  Cumulative  worldwide  installed  photovoltaics  reaches  1000  MW 

•  20(H)  Olympics  in  Australia  highlight  wide  range  of  PV  applications,  and  the  awarding  of  the  lirsl 
Bachelor  of  Engineering  degrees  in  Photovoltaics  and  Solar  Engineering  (UNSW.  Australia) 

•  2(M)2  Cumulative  worldwide  installed  photovoltaics  reaches  2(HK)  MW.  It  took  25  years  to  reach  the 
lirsl  l(HM)  MW  and  only  3  years  to  double  it:  production  of  crystalline  Si  cells  exceeds  100  MW  per 
year  at  Sharp  Corp.  (Japan).  BP  Solar  ceases  RAD  and  production  of  a-Si  and  CdTe  thin-lilm  modules 
in  USA  ending  >20  years  of  effort 

Table  3.  Notable  events  in  the  history  of  photovoltaics  (From 

[12]). 
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The  Electromagnetic  Spectrum 


a . 


Light  is  electromagnetic  radiation.  The  frequency 
of  light  determines  its  color.  Figure  11  shows  the  visible 
part  of  the  electromagnetic  spectrum.  Visible  wavelengths 
range  from  390  nm  (violet)  to  780  nm  (red) . 


VISIBLE 


Wavelength  (m) 

K>«  10” 

r«  105  10* 

1  1  m 

Ultraviolet 

1  1 

1  1 

Infrared 

1  l 

10' 7  10*  K 

Frequency  (Hz) 

i  i  l 

*5 

1i 

w  to15  10” 

KP  10J  10' 

Photon  energy  (eV) 

io> 

lO'  10” 

SOLAR 

RADIATION 


Figure  11 


The  electromagnetic  spectrum  (From  [13]) 


Table  4  shows  the  approximate  wavelength  range  of 
visible  colors. 


Color 

Wavelength  (nm) 

Red 

780-622 

Orange 

622-597 

Yellow 

597-577 

Green 

577-492 

Blue 

492-455 

Violet 

455-390 

Table  4 .  Approximate  wavelength  of  various  colors  in  vacuum 

(After  [  14 ] )  . 
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The  sun  emits  light  from  ultraviolet,  visible,  and 
infrared  wavelengths  in  the  electromagnetic  spectrum.  Solar 
irradiance  has  the  largest  magnitude  at  visible  wavelengths, 
peaking  in  the  blue-green  [10,  p.  17] . 

b .  Band  Gap 

The  band  gap  of  the  semiconductor  material 
determines  how  the  solar  cell  reacts  to  light.  Table  5  shows 
a  small  sample  of  semiconductor  band  gaps.  Chapter  IV  covers 
the  band  gaps  of  Indium  Nitride,  Gallium  Nitride,  and  Indium 
Gallium  Nitride. 


Material 

Band  gap  (eV)  at  300  K 

Si 

1 . 12 

Ge 

0.66 

GsAs 

1 . 42 

InP 

1.34 

Table  5.  Common  semiconductor  band  gaps  (After  [9,  p.  31]) . 


The  band  gap  of  the  semiconductor  material 
determines  the  wavelength  of  light  that  meet  the 
requirements  to  generate  electrical  energy.  The  conversion 
formula  between  band  gap  and  wavelength  is: 

A(jum)  =  — — — 

Eg(eV) 

1  24 

A  (jum)  = - 

Eg(eV) 

Where  A  is  the  wavelength  in  micrometers,  h  is 
Planck' s  constant,  c  is  the  speed  of  light  in  vacuum,  and  Eg 
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is  the  band  gap  in  eV.  One  eV  is  approximately  equal  to 
1.6xl0“19  J  of  energy.  In  the  case  of  Gallium  Arsenide,  the 
wavelength  that  corresponds  to  1.42  eV  is  0.873  fim  . 

Figure  12  helps  visualize  the  concept  of  light 
absorption.  When  light  has  energy  greater  than  1.1  eV,  the 
silicon  solar  cell  generates  electricity.  Light  with  less 
than  1.1  eV  of  energy  is  unused.  Similarly,  light  with 
energy  greater  than  1.43  eV  excites  the  outer  shell 
electrons  of  the  gallium  arsenide  solar  cell.  And  finally, 
light  with  energy  greater  than  1.7  eV  is  useful  for  aluminum 
gallium  arsenide  photovoltaic  material. 


Silicon  Gallium  arsenide  Aluminum 

absorbs  absorbs  gallium  arsenide 

>  l.l  eV  >l41eV  absorbs 

>1.7  eV 


<  l.l  eV  <  1 .43  eV  <l.7eV 


Figure  12.  Effect  of  light  energy  on  different  band  gaps 

(From  [  15 ] )  . 


c.  Solar  Cell  Junctions 

The  discussion  in  the  previous  section  discussed 
the  effect  of  light  energy  on  different  band  gaps.  Treated 
individually,  each  of  the  photovoltaic  materials  from  Figure 
12  would  act  as  a  single  junction  solar  cell. 

However,  to  increase  the  efficiency  of  the  solar 
cell,  multiple  junctions  can  be  created.  For  example,  in 
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Figure  12,  the  top  junction  would  be  made  up  of  Aluminum 
Gallium  Arsenide.  This  junction  would  absorb  light  energy 
greater  than  1.7  eV.  Any  unused  photons  would  be  filtered 
through  to  the  next  junction.  The  gallium  arsenide  junction 
would  then  absorb  the  photons  with  energy  greater  than  1.4 
eV.  The  remaining  photons  would  be  absorbed  by  the  silicon 
junction . 

Although  the  above  paragraph  described  the  basics 
of  a  multi j unction  solar  cell,  such  device  may  not  produce 
the  desired  results  due  to  lattice  mismatch.  The  next 
section  covers  the  basics  of  lattice  matching. 

d.  Lattice  Matching 

Semiconductors  are  three-dimensional  in  their  cell 
structure.  The  simple  cubic  structure  serves  to  illustrate 
the  concept  of  lattice  and  is  presented  in  Figure  13. 


Figure  13.  Simple  cubic  lattice  structure  (From  [16]) . 


Figure  13  shows  that  each  side  of  the  cube  is 
represented  by  the  letter  "a".  The  separation  "a"  is  called 
the  lattice  constant.  Each  semiconductor  material  has  a 
lattice  constant.  Therefore,  when  creating  multi j unction 
solar  cells,  the  lattices  must  be  matched. 
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Figure  14  shows  the  lattice  constants  for  several 
semiconductors.  An  example  given  by  P.  Michalopoulos  [1,  p. 
87]  shows  how  to  lattice  match  Indium  Gallium  Phosphide  to 
Gallium  Arsenide. 


The  alloy  Indium  Gallium  Phosphide  is  composed  of 
x  parts  of  Gallium  phosphide  and  1-x  parts  of  Indium 
Phosphide.  Therefore,  Indium  Gallium  Phosphide  is 
represented  as  Ini_xGaxP. 

GaAs  has  a  lattice  constant  a=5.65A,  GaP  has 
a=5.45A  and  InP  has  a=5.87A.  The  goal  is  to  create  InGaP 
with  a  lattice  constant  that  matches  that  of  GaAs.  The 
formula  to  find  x  is  given  as: 
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aGaAs  =  aGaP  '  *  +  ainP  ’  (1 "  *)  <=>  X 


_  alnP 

aGaP  —  ainP 

With  x=0.52,  Ino.48Gao.52P  has  a=5.65A.  A  rough 
approximation  of  the  resulting  band  gap  of  Ino.4sGao.52P  is 
given  as  follows: 

E^GaP=xEgaP+(l-x)Egp 

The  eguation  yields  a  band  gap  of  1.9  eV  for 
Ino.4sGao.52P •  Therefore,  a  dual- j unction  solar  cell  of  InGaP 
at  1.9  eV  and  GaAs  at  1.4  eV  can  be  built. 

Figure  15  shows  the  lattice  constant  for  Indium 
Nitride  and  Gallium  Nitride. 


Figure  15.  Lattice  constant  for  InN  and  GaN  (From  [18])  . 

e .  AMO  Spectrum 

The  location  of  the  solar  cell  affects  the  input 
solar  radiation  spectrum.  A  solar  cell  on  Mars  receives  a 
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different  (smaller)  spectrum  than  a  solar  cell  on  a 
satellite  that  orbits  Earth.  The  energy  received  outside 
Earth's  atmosphere  is  approximately  1365  W/m2.  This  spectrum 
is  called  Air  Mass  Zero  or  AMO.  Terrestrial  solar  cells  have 
to  deal  with  the  attenuation  of  the  solar  spectrum  due  to 
Earth's  atmosphere.  This  solar  spectrum  is  called  AMI . 5 .  For 
the  purposes  of  this  thesis,  AMO  is  used  during  simulations. 
From  data  obtained  from  the  National  Renewable  Energy 
Laboratory  (NREL)  [19],  the  AMO  spectrum  was  plotted  using  a 
Matlab  script. 

Figures  16  and  17  show  the  AMO  spectrum  with 
respect  to  wavelength  and  energy,  respectively.  From  Figure 
17,  it  can  be  seen  that  semiconductor  materials  with  band 
gaps  of  less  than  4  eV  are  able  to  extract  most  of  the  solar 
spectrum. 


AMO 


Figure  16.  AMO  spectrum  (Wavelength  vs  Irradiance)  (After 

[19]). 
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AMO 


Figure  17.  AMO  spectrum  (Energy  vs  Irradiance)  (After  [19]) . 

f.  Current-Voltage  Curves 

A  typical  solar  cell  current-voltage  (IV)  curve  is 
presented  in  Figure  18. 


Cell  voltage 

[V] 

Figure  18.  Sample  IV  curve  used  in  efficiency  calculations 

(After  [20] ) . 
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From  Figure  18,  there  are  several  points  of 
interest.  The  short  circuit  current  (ISc)  occurs  when  the 
voltage  is  zero.  This  is  the  highest  absolute  value  current. 
The  open  circuit  voltage  (Voc)  occurs  when  the  current  is 
zero.  This  is  the  highest  voltage.  The  dimensions  of  the 
larger  rectangle  in  Figure  18  are  determined  by  V0c  and  Isc. 
Since  power  (P)  is  determined  by  the  product  of  current 
times  voltage,  the  maximum  power  point  occurs  at  (VMP,  Imp)  • 

The  calculations  for  solar  cell  efficiency  are  as 

follows : 


P  max — Imp  V  mp 

P  max  Imp V mp 


FF= 
rt 


Isc  V  oc  Isc  V  oc 

_  P  max _ Imp  V mp 

=  ~PhT_  Pin 


Where  Pmax  is  the  maximum  power  point,  FF  is  the 
fill  factor,  and  77  is  the  efficiency.  The  fill  factor 
measures  the  "squareness"  of  the  IV  curve. 


g.  Electrical  Output 

A  solar  cell  is  a  p-n  junction  photodiode.  In 
order  to  obtain  the  IV  characteristic  of  the  solar  cell,  the 
dark  current  needs  to  be  subtracted  from  the  photogenerated 
current . 

1  =  II  — Id 

The  dark  current  is  the  current  through  the  solar 
cell  when  bias  is  applied  in  the  dark  [10,  p.  30] . 

Graphically,  the  IV  characteristic  is  obtained  in 

Figure  19. 
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Figure  19.  Solar  cell  IV  characteristic  (From  [21]). 


C.  CHAPTER  CONCLUSIONS 

This  chapter  provided  basic  information  on 
semiconductors  and  solar  cells.  The  foundation  has  been 
established  from  the  physics  stand  point.  The  next  step  is 
to  cover  an  introduction  to  the  simulation  software. 
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Ill .  SILVACO  ATLAS  SIMULATION  SOFTWARE 


This  thesis  uses  Silvaco  Atlas  to  perform  solar  cell 
modeling.  Silvaco  International  produces  a  suite  of  software 
programs  that  together  become  a  Virtual  Wafer  Fabrication 
tool.  This  chapter  introduces  Silvaco  Atlas  and  some  of  its 
features . 


A.  VIRTUAL  WAFER  FAB 


Silvaco  International  provides  several  software  tools 
to  perform  process  and  device  simulation.  From  Figure  20,  it 
can  be  seen  that  Silvaco  offers  powerful  simulation 
software . 


ATHENA  I  SSuprem3 

1D/2D  Process  I  ID  Process 
Simulation  Software  I  Simulation  Software 


ATLAS  MERCURY 

Device  Simulation  FastFET  Simulation 
Software  Software 


III 


|  CLEVER 

Monts  Cstlo  I  Physics-bassd 

Transport  Parameter  I  Parasitic  Extraction 
Generator  I  Software 


it  l 


UTMOST 

SPICE  Modeling 


Virtual  Wafer  Fab  Interactive  Environment 


DeckBuild  I  DevEdit  I  TonyPlot  1D/2D  I  TonyPlot3D  I  C-Interpreter  I  n  ti  •  I  MackViPwc 

Run-Time  I  Structure  &  Mesh  I  Interactive  Visualization  I  Interactive  Visualization  I  User  Defined  I  utility  I  Lavout  Editor 
Environment  I  Editor  I  Utility  I  Utility  I  Models  I  ty  ■  y 


Platforms:  Solaris,  HPUX,  Linux  and  Windows 


Figure  20.  Silvaco' s  Virtual  Wafer  Fabrication  Environment 

(From  [22 ] )  . 


In  this  thesis,  Silvaco  Atlas  was  extensively  used.  The 
DeckBuild  run-time  environment  received  the  input  files. 
Within  the  input  files,  Silvaco  Atlas  was  called  to  execute 
the  code.  And  finally,  TonyPlot  was  used  to  view  the  output 
of  the  simulation.  Additionally,  output  log  files  were 
produced.  The  data  extracted  from  the  log  files  could  then 
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be  displayed  using  Microsoft  Excel  or  Matlab  scripts.  Figure 
21  shows  the  inputs  and  outputs  for  Silvaco  Atlas. 


Figure  21.  Atlas  inputs  and  outputs  (From  [23,  p.  2-2]). 

B.  SILVACO  ATLAS 

Atlas  is  a  software  program  used  to  simulate  two  and 
three-dimensional  semiconductor  devices.  The  physical  models 
include  in  Atlas  are  presented  in  Table  6. 
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•  DC,  AC  small-signal,  and  full  time-dependency. 

•  Drift -diffusion  transport  models. 

•  Energy  balance  and  Hydrodynamic  transport  models. 

•  Lattice  heating  and  heatsinks. 

•  Graded  and  abrupt  heterojunctions. 

•  Optoelectronic  interactions  with  general  ray  tracing. 

•  Amorphous  and  polycrvstalline  materials. 

•  General  circuit  environments. 

•  Stimulated  emission  and  radiation 

•  Fermi-Dirac  and  Boltzmann  statistics. 

•  Advanced  mobility  models. 

•  Heavy  doping  effects. 

•  F  ull  acceptor  and  donor  trap  dynamics 

•  Ohmic,  Schottky,  and  insulating  contacts. 

•  SRH,  radiative,  Auger,  and  surface  recombination. 

•  Impact  ionization  (local  and  non-local). 

•  Floating  gates. 

•  Band-to-band  and  Fowler-Nordheim  tunneling. 

•  Hot  carrier  injection. 

•  Quantum  transport  models 

•  Thermionic  emission  currents. 

Table  6.  Silvaco  Atlas  physical  models  (From  [23,  p.  1-2]). 


C.  INPUT  FILE  STRUCTURE 

Silvaco  Atlas  receives  input  files  through  DeckBuild. 
The  code  entered  in  the  input  file  calls  Atlas  to  run  with 
the  following  command: 

go  atlas 

Following  that  command,  the  input  file  needs  to  follow 
a  pattern.  The  command  groups  are  listed  in  Figure  22. 
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Group 


Statements 


1.  Structure  Specification 


MESH 

REGION 

ELECTRODE 

DOPING 


2.  Material  Models  Specification 


MATERIAL 

MODELS 

CONTACT 

INTERFACE 


3.  Numerical  Method  Selection  METHOD 


4.  Solution  Specification 


LOG 

SOLVE 

LOAD 

SAVE 


5.  Results  Analysis  EXTRACT 

TONYPLOT 


Figure  22.  Atlas  command  groups  and  primary  statements  (From 

[23,  p.  2-8] ) . 


Atlas  follows  the  following  format  for  statements  and 
parameters : 

<STATEMENT>  <PARAMETER>=<VALUE> 

The  following  line  of  code  serves  as  an  example. 

DOPING  UNIFORM  N . TYPE  CONCENTRATIONS . 0el6  REGIONS  \ 
OUTFILE=my . dop 

The  statement  is  DOPING.  The  parameters  are  UNIFORM, 
N . TYPE ,  CONCENTRATION,  REGION,  and  OUTFILE.  There  are  four 
different  type  of  parameters:  real,  integer,  character,  and 
logical.  The  back  slash  (\)  serves  the  purpose  of  continuing 
the  code  in  the  next  line.  Parameters,  such  as  UNIFORM,  are 
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logical.  Unless  a  TRUE  or  FALSE  value  is  assigned,  the 

parameter  is  assigned  the  default  value.  This  value  can  be 
either  TRUE  or  FALSE.  The  Silvaco  Atlas  manual  needs  to  be 
referenced  to  identify  the  default  value  assigned  to 
specific  parameters. 

1.  Structure  Specification 

The  structure  specification  is  done  by  defining  the 

mesh,  the  region,  the  electrodes  and  the  doping  levels. 

a .  Mesh 

The  mesh  used  for  this  thesis  is  two-dimensional. 
Therefore,  only  x  and  y  parameters  are  defined.  The  mesh  is 
a  series  of  horizontal  and  vertical  lines  and  spacing 
between  them.  From  Figure  23,  the  mesh  statements  are 

specified. 


-2.8 


-2.4 


-2 


-1.6 


-1.2 


-0.8 
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-200 


“1 — 

-100 


0 

Microns 


x.mesh 
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x.mesh 

loc=-25 
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x.mesh 
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x.mesh 

loc=2S0 

spac=25 

y.mesh 

loc=-2. 9 

spac=0. 01 

y.mesh 

loc=-2. 8 

spac=0. 03 

y.mesh 

loc=-2. 5 

spac=0. 05 

y.mesh 

loc=-2 

spac=0. 05 

y.mesh 

loc=-l 

spac=0. 1 

y.mesh 
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spac=0. 1 

~r~ 

100 


200 


Figure  23 


Atlas  mesh  (From  2,  p.18]) 
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The  general  format  to  define  the  mesh  is: 


X .  MESH  LOCATION=<VALUE>  SPACING=<VALUE> 

Y .  MESH  LOCATION=<VALUE>  SPACING=<VALUE> 

For  example,  the  x.mesh  starting  at  -250  microns 
has  spacing  of  25  microns.  That  means  it  is  relatively 
coarse.  The  x.mesh  becomes  finer  between  -25  and  25  microns 
with  a  spacing  of  2.5  microns.  The  y.mesh  is  similarly 
defined.  For  example,  at  y.mesh  of  -2.9  microns,  the  spacing 
is  0.01  microns.  Then  at  location  y.mesh  of  -2.8  microns, 
the  spacing  changes  to  0.03  microns.  The  mesh  is  coarser  at 
y.mesh  location  of  -1,  when  the  spacing  is  0.1. 

A  coarse  or  fine  mesh  determines  the  accuracy  of 
the  simulation.  A  coarse  mesh  produces  a  faster  simulation, 
but  less  accurate  results.  A  fine  mesh  produces  a  slower 
simulation,  but  more  accurate  results.  The  areas  that  have  a 
finer  mesh,  therefore,  are  of  greatest  interest  in  the 
simulation . 


b.  Region 

After  defining  the  mesh,  it  is  necessary  to  define 
the  regions.  The  format  to  define  the  regions  is  as  follows: 

REGION  number=<integer>  <material_type>  / 

<position  parameters> 

From  Figure  24,  the  code  that  defines  the  regions 
is  identified.  There  are  six  regions  defined.  The  limits  of 
each  region  are  explicitly  identified  in  the  x-  and  y-axis. 
The  regions  must  then  be  given  a  material. 
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Figure  24.  Atlas  region  (From  [2,  p.  19]) 


From  Figure  25,  the  code  defines  the  material  for 
each  region.  Note  that  the  color  coding  identifies  the 
material.  The  regions  have  vertical  and  horizontal  lines  to 
mark  their  boundaries. 


-2.4 


S-1.6 


-1.2 


-0.0 


-0.4 


Figure  25.  Atlas  regions  with  materials  defined  (From  [2,  p. 

19]). 
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c.  Electrodes 

The  next  structure  specification  corresponds  to 
electrodes.  Typically,  in  this  simulation  the  only 
electrodes  defined  are  the  anode  and  the  cathode.  However, 
Silvaco  Atlas  has  a  limit  of  50  electrodes  that  can  be 
defined.  The  format  to  define  electrodes  is  as  follows: 

ELECTRODE  NAME=<electrode  name>  <position_parameters> 

From  Figure  26,  the  electrode  statements  are 
defined  for  the  anode  and  the  cathode.  Note  that  the  cathode 
is  defined  with  gold  as  the  material.  The  x  and  y  dimensions 
correspond  to  region  6  previously  defined.  Meanwhile,  the 
anode  is  defined  at  the  bottom  of  the  cell  for  the  entire  x- 
range  at  y=0 . 


electrode  name=cathode  material =Gold  x.min=-25  x.max=25  y.min=-2.9  y.max=-2.8 
electrode  name=anode  x.min=-250  x.max=250  y.min=0  y.max=0 


-200  -100  0  100  200 
Microns 


Figure  26.  Atlas  electrodes  (From  [2,  p.  20]). 
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d.  Doping 


The  last  aspect  of  structure  specification  that 
needs  to  be  defined  is  doping.  The  format  of  the  Atlas 
statement  is  as  follows: 

DOPING  distribution  type>  <dopant_type>  / 

<position  parameters> 

From  Figure  27,  the  doping  types  and  the  doping 
levels  are  defined.  Doping  can  be  n-type  or  p-type.  The 
distribution  type  can  be  uniform  or  Gaussian. 


-0.8 


-0.4 


doping  uniform  region=l  n. type  conc=2el8 
doping  uniform  region=4  n. type  conc=le20 
doping  uniform  region=5  p.type  conc=8e!8 


n — I — r 

o 

Microns 


I — 
100 


I — 

200 


Figure  27.  Atlas  doping  (From  [2,  p.  21]). 


2 .  Materials  Model  Specification 

After  the  structure  specification,  the  materials  model 
specification  is  next.  From  Figure  28,  the  materials  model 
specification  is  broken  down  into  material,  models,  contact, 
and  interface. 
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2.  Material  Models  Specification 


MATERIAL 
MODELS 
CONTACT 
INTERFACE 

Figure  28.  Atlas  material  models  specification  (After  [23,  p. 

2-8]). 

a.  Material 

The  format  for  the  material  statement  is  as 

follows : 

MATERIAL  <localization>  <material_def inition> 

Below  are  three  examples  of  the  material 

statement : 

MATERIAL  MATERIAL=Silicon  EG300=1.1  MUN=1200 

MATERIAL  REGI0N=4  TAUN0=3e-7  TAUP0=2e-5 

MATERIAL  NAME=base  NC300=4el8 

In  all  examples,  when  MATERIAL  appears  first,  it 
is  considered  the  statement.  When  MATERIAL  appears  a  second 
time  in  the  first  example,  it  is  considered  a  localization 
parameter.  In  the  second  and  third  examples,  the 
localization  parameters  are  REGION  and  NAME,  respectively. 
Various  other  parameters  can  be  defined  with  the  material 
statement.  Examples  of  these  parameters  are  the  band  gap  at 
room  temperature  (EG300) ,  electron  mobility  (MUN) ,  electron 
(TAUNO)  and  hole  (TAUPO)  recombination  lifetimes,  conduction 
band  density  at  room  temperature  (NC300) ,  among  others. 

Jb.  Models 

The  physical  models  fall  into  five  categories: 
mobility,  recombination,  carrier  statistics,  impact 
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ionization,  and  tunneling.  The  syntax  of  the  model  statement 
is  as  follows: 

MODELS  <model  flag>  <general  parameter>  / 

<model  dependent  parameters> 

The  choice  of  model  depends  on  the  materials 
chosen  for  simulation. 

The  example  below  activates  several  models. 

MODELS  CONMOB  FLDMOB  SRH 

CONMOB  is  the  concentration  dependent  model. 
FLDMOB  is  the  parallel  electric  field  dependence  model.  SRH 
is  the  Shockley-Read-Hall  model. 

c.  Contact 

Contact  determines  the  attributes  of  the 
electrode.  The  syntax  for  contact  is  as  follows: 

CONTACT  NUMBER=<n>  | NAME=<ename> | ALL 

The  following  is  an  example  of  the  contact 

statement . 

CONTACT  NAME=anode  current 

d.  Interface 

The  semiconductor  or  insulator  boundaries  are 
determined  with  the  interface  statement.  The  syntax  is  as 
follows : 

INTERFACE  [ <parameters> ] 

The  following  example  shows  the  usage  of  the 
interface  statement. 
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INTERFACE  X.MIN=-4  X.MAX=4  Y.MIN=-0.5  Y.MAX=4  \ 


QF=lel 0  S . N=le4  S.P=le4 

The  max  and  min  values  determine  the  boundaries. 
The  QF  value  specifies  the  fixed  oxide  charge  density  (cm' 
2)  .  The  S.N  value  specifies  the  electron  surface 
recombination  velocity.  S.P  is  similar  to  S.N,  but  for 
holes . 


3 .  Numerical  Method  Selection 

After  the  materials  model  specification,  the  numerical 
method  selection  must  be  specified.  From  Figure  29,  the  only 
statement  that  applies  to  numerical  method  selection  is 
method . 

3.  Numerical  Method  Selection  METHOD 

Figure  29.  Atlas  numerical  method  selection  (After  [23,  p.  2- 

8  ]  )  . 

There  are  various  numerical  methods  to  calculate 
solutions  to  semiconductor  device  problems.  There  are  three 
types  of  solution  techniques  used  in  Silvaco  Atlas: 

•  decoupled  (GUMMEL) 

•  fully  coupled  (NEWTON) 

•  BLOCK 

The  GUMMEL  method  solves  for  each  unknowns  by  keeping 
all  other  unknowns  constant.  The  process  is  repeated  until 
there  is  a  stable  solution.  The  NEWTON  method  solves  all 
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unknowns  simultaneously.  The  BLOCK  method  solves  some 
equations  with  the  GUMMEL  method  and  some  with  the  NEWTON 
method . 

The  GUMMEL  method  is  used  for  a  system  of  equations 
that  are  weakly  coupled  and  there  is  linear  convergence.  The 
NEWTON  method  is  used  when  equations  are  strongly  coupled 
and  there  is  quadratic  convergence. 

The  following  example  shows  the  use  of  the  method 
statement . 

METHOD  GUMMEL  NEWTON 

In  this  example,  the  equations  are  solved  with  the 
GUMMEL  method.  If  convergence  is  not  achieved,  then  the 
equations  are  solved  using  the  NEWTON  method. 


4.  Solution  Specification 


After  completing  the  numerical  method  selection,  the 
solution  specification  is  next.  Solution  specification  is 
broken  down  into  log,  solve,  load,  and  save  statements,  as 
shown  in  Figure  30. 


4.  Solution  Specification 


LOG 

SOLVE 

LOAD 

SAVE 


Figure  30.  Atlas  solution  specification  (After  [23,  p.  2-8]). 


a .  Log 


LOG  saves  all  terminal  characteristics  to  a  file. 
DC,  transient,  or  AC  data  generated  by  a  SOLVE  statement 
after  a  LOG  statement  is  saved. 
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The  following  shows  an  example  of  the  LOG 

statement . 

LOG  OUTFILE=myoutputf ile . log 

The  example  saves  the  current-voltage  information 
into  myoutputfile.log. 

b.  Solve 

The  SOLVE  statement  follows  the  LOG  statement. 
SOLVE  performs  a  solution  for  one  or  more  bias  points.  The 
following  is  an  example  of  the  SOLVE  statement. 

SOLVE  Bl=l 0  B3=5  BEAM=1  SS.PHOT  SS . LIGHT=0 . 0 1  \ 

MULT . F  FREQUENCY=le3  FSTEP=10  NFSTEP=6 


B1  and 

B3 

specify 

the 

optical 

spot 

power 

associated  with 

the 

optical 

beam 

numbers 

1  and  3, 

respectively.  The 

beam 

number  is 

an 

integer  between 

1  and 

10.  BEAM  is  the  beam  number  of  the  optical  beam  during  AC 
photogeneration  analysis.  SS.PHOT  is  the  small  signal  AC 
analysis.  SS. LIGHT  is  the  intensity  of  the  small  signal  part 
of  the  optical  beam  during  signal  AC  photogeneration 
analysis.  MULT . F  is  the  frequency  to  be  multiplied  by  FSTEP. 
NFSTEPS  is  the  number  of  times  that  the  frequency  is 
incremented  by  FSTEP. 

c.  Load  and  Save 

The  LOAD  statement  enters  previous  solutions  from 
files  as  initial  guess  to  other  bias  points.  The  SAVE 
statement  enters  all  node  point  information  into  an  output 
f  ile . 
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The  following  are  examples  of  LOAD  and  SAVE 
statements . 

SAVE  OUTF=SOL . STR 

In  this  case,  SOL. STR  has  information  saved  after 
a  SOLVE  statement.  Then,  in  a  different  simulation,  SOL. STR 
can  be  loaded  as  follows: 

LOAD  INFILE=SOL . STR 

5 .  Results  Analysis 

Once  a  solution  has  been  found  for  a  semiconductor 
device  problem,  the  information  can  be  displayed  graphically 
with  TonyPlot.  Additionally,  device  parameters  can  be 
extracted  with  the  EXTRACT  statement,  as  shown  in  Figure  31. 

5.  Results  Analysis  EXTRACT 

TONYPLOT 

Figure  31.  Atlas  results  analysis  (After  [23,  p.  2-8]). 

In  the  example  below,  the  EXTRACT  statement  obtains  the 
current  and  voltage  characteristics  of  a  solar  cell.  This 
information  is  saved  into  the  IVcurve.dat  file.  Then, 
TonyPlot  plots  the  information  in  the  IVcurve.dat  file. 

EXTRACT  NAME="iv"  curve (v. "anode" ,  I. "cathode")  / 

OUTFILE=" IVcurve . dat " 

TONYPLOT  IVcurve.dat 

Figure  32  shows  the  sample  IV  curve  plotted  by 
TonyPlot . 
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Figure  32.  Sample  TonyPlot  IV  curve. 


D .  CONCLUSION 


This  chapter  presented  an  introduction  to 
Atlas,  the  structure  of  the  input  files,  and  some 
statements.  With  these  basic  tools,  the  next 
introduces  some  basic  information  on  wurtzite  Indium 
Nitride  before  proceeding  to  the  thesis  simulation. 


Silvaco 
of  its 
chapter 
Gallium 
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IV.  INDIUM  GALLIUM  NITRIDE 


A.  A  FULL  SPECTRUM  PHOTOVOLTAIC  MATERIAL 

Prior  to  2001,  it  was  thought  that  wurtzite  Indium 
Nitride  had  a  band  gap  of  approximately  1.9  eV.  In  2001,  it 
was  discovered  that  Indium  Nitride  had  a  much  smaller  band 
gap.  Davydov,  et  al . ,  concluded  in  [24]  that  Indium  Nitride 
had  a  band  gap  of  approximately  0.9  eV.  Additionally, 
Davydov,  et  al .  ,  presented  in  [25]  that  wurtzite  Indium 
Gallium  Nitride  band  gaps  in  the  range  x=0.36  to  x=l 
supported  findings  in  [24] .  The  band  gaps  for  Indium  Gallium 
Nitride,  with  Indium  Nitride  concentrations  ranging  from  0.0 
to  1.0  are  presented  in  Figure  33. 


InN  fraction,  x 

Figure  33.  InGaN  band  gap  as  a  function  of  In  composition 

(After  [25] )  . 
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Later  on,  research  conducted  at  Lawrence  Berkeley 
National  Laboratory  presented  that  the  Indium  Nitride 
fundamental  band  gap  was  approximately  0.77  eV  at  room 
temperature.  Since  Gallium  Nitride  has  a  band  gap  of 
approximately  3.4  eV  at  room  temperature,  then  Indium 
Gallium  Nitride  can  have  a  band  gap  ranging  from  0.77  eV  to 
3.4  eV  by  changing  the  percent  composition  of  Indium  and 
Gallium  within  Indium  Gallium  Nitride.  Figure  34  confirms 
that  Indium  Gallium  Nitride  follows  the  pattern  of  ranging 
from  0.77  eV  to  3.4  eV. 


(1021  photons/sec/mVn)  ln11t  Ga^N  alloy 

Figure  34.  InGaN  band  gap  and  solar  spectrum  comparison 

(After  [26]) . 

Figure  34  also  presents  the  solar  spectrum  on  the  left- 
hand  side.  Therefore,  it  is  clear  that  the  visible  spectrum 
has  a  near-perfect  match  with  the  Indium  Gallium  Nitride 
range  of  band  gaps . 

It  should  be  noted  that  research  conducted  in  [27] 
presents  a  Indium  Nitride  band  gap  of  1.7 +  0.2  eV.  It  is 
stated  in  [27],  "These  investigations  clarify  that  a  band 
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transition  around  1.7  eV  does  exist  in  wurtzite  InN  although 
in  this  energy  range  optical  material  responses  are  hardly 
ever  found.  If  the  1.7  eV  band  transition  is  indeed  the 
fundamental  bandgap  in  InN  at  room  temperature  it  is 
concluded  that  defect  bands,  grain  boundaries,  dislocations 
and/or  the  very  conductive  surface  are  contributing  to  the 
lower  energy  optical  responses  in  InN.  The  identification  of 
the  various  defects  in  InN  will  be  subject  of  future 
investigations . " 

In  follow  up  research,  [28]  explains  the  possible 
reason  for  the  higher  Indium  Nitride  band  gap:  "Owing  to 
InN' s  exceptional  propensity  for  n-type  activity,  in  many 
instances  the  larger  apparent  values  of  the  band  gap  can  be 
attributed  to  the  Burstein-Moss  shift  of  the  optical 
absorption  edge  resulting  from  the  occupation  of  conduction 
band  states  by  free  electrons." 

Figure  35  shows  evidence  that  the  fundamental  band  gap 
of  indium  nitride  is  approximately  0.7  eV,  and  not  near  2 
eV.  The  optical  absorption  (blue  curve)  has  an  onset  at 
approximately  0.8  eV.  There  is  no  change  in  absorption 
around  the  2  eV  range.  The  photo  luminescence  (red  curve) 
shows  a  peak  at  the  band  edge  of  approximately  0.8  eV. 
Finally,  the  photo-modulated  reflectance  (green  curve) 
presents  a  transition  close  to  0.8  eV.  These  three  curves 
support  a  band  gap  for  indium  nitride  of  approximately  0.8 
eV.  The  1 . 9-2 . 0  eV  band  gap  is  not  evident  in  Figure  35. 
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Figure  35. 


Evidence  of  0.7  eV  band  gap  for  indium  nitride 
(From  [29]) . 


The  controversy  about  indium  nitride's  fundamental  band 
gap  is  not  the  subject  of  this  thesis.  This  research 
proceeded  accepting  the  fundamental  band  gap  of  0.77  eV  for 
indium  nitride. 

According  to  [30],  the  following  formula  provides  an 
approximation  of  Indium  Gallium  Nitride  band  gap: 

Eg(x)  =  3 .42x  +  0.77(1  -  x)  - 1 .43x(l  -  x) 

Where  EG  is  the  InGaN  band  gap,  3.42  eV  is  the  GaN  band 
gap,  0.77  eV  is  the  InN  band  gap,  1.43  eV  is  the  bowing 
parameter  b,  x  is  the  Ga  concentration,  and  (1-x)  is  the  In 
concentration . 

When  plotting  the  band  gap  formula  for  InGaN  using  a 
Matlab  script,  the  Figure  36  is  obtained. 
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Indium  concentration  vs  InGaN  Band  Gap 


Figure  36.  InGaN  band  gap  as  a  function  of  In  concentration. 


The  band  gaps  used  in  the  Silvaco  Atlas  simulations 
were  found  by  modifying  the  Matlab  script  for  Figure  35. 

B.  RADIATION-HARD  SEMICONDUCTOR  MATERIAL 

Radiation-hard  materials  are  essential  for  space 
applications.  When  compared  to  other  Gallium  Arsenide  and 
Indium  Gallium  Phosphide,  Indium  Gallium  Nitride  is  able  to 
withstand  a  greater  amount  of  radiation.  According  to  [31], 
InGaN  "retains  its  optoelectronic  properties  at  radiation 
damage  doses  at  least  two  orders  of  magnitude  higher  than 
the  damage  thresholds  of  the  materials  (GaAs  and  GalnP) 
currently  used  in  high  efficiency  MJ  cells".  Radiation 
levels  in  [31]  were  1  MeV  electron,  2  MeV  proton,  and  2  MeV 
alpha  particle  irradiation. 
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This  indicates  that  InGaN  is  not  only  potentially  a 
high-efficiency  photovoltaic  material,  but  it  is  also  able 
to  withstand  the  harsh  space  environment. 

C.  INDIUM  GALLIUM  NITRIDE  CHALLENGES 

There  are  some  significant  issues  that  currently 
prevent  InGaN  from  being  used  as  photovoltaic  material.  As 
of  this  writing,  no  PN  junction  has  been  built.  The  PN 
junction  is  essential  to  create  photodiodes.  Indium  Nitride 
has  the  highest  electron  affinity  among  all  semiconductors 
[32] .  Hence,  Indium  Nitride  has  a  tendency  to  be  n-type.  It 
is  much  more  difficult  to  build  p-type  material  with  Indium 
Nitride.  However,  evidence  of  p-type  doping  in  Indium 
Nitride  and  Indium  Gallium  Nitride  has  been  reported  in  [33] 
using  Magnesium.  This  is  a  step  in  the  right  direction. 
Further  progress  in  this  area  dictates  the  creation  of 
Indium  Gallium  Nitride  PN  junctions. 

Having  reviewed  the  characteristics  of  Indium  Gallium 
Nitride,  the  next  chapter  covers  the  simulation  of  Indium 
Gallium  Nitride  using  Silvaco  Atlas. 
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V.  SIMULATION  OF  INDIUM  GALLIUM  NITRIDE  IN  SILVACO 

ATLAS 


One  of  the  primary  motivations  behind  this  research  was 
to  find  a  way  to  make  solar  cells  significantly  more 
efficient.  When  researching  high-efficiency  solar  cells,  the 
improvements  were  typically  incremental  over  previous 
results.  Solar  cells  were  first  designed  as  single- j unction . 
The  most  popular  photovoltaic  material  has  been  Silicon  due 
to  its  abundance  and  relatively  low  cost.  The  maximum 
efficiency  of  polycrystalline  Silicon  is  approximately 
19.8%.  Commercial  silicon  solar  cells  have  a  maximum 
efficiency  of  approximately  13%.  Improvements  in  efficiency 
can  be  obtained  in  single- j unction  solar  cells  by  using 
Gallium  Arsenide.  Monocrystalline  Gallium  Arsenide  has  a 
maximum  efficiency  of  approximately  25.1%  [10,  p.  195]. 

By  designing  dual- j unction  solar  cells,  the  efficiency 
can  again  be  improved.  A  dual- j unction  cell  with  Gallium 
Phosphide  as  the  top  junction  and  Gallium  Arsenide  as  the 
bottom  cell  can  have  an  efficiency  of  30.3%  [10,  p  302] . 

Increasing  the  number  of  junctions  makes  the  cells  more 
complex  in  design  and  the  gains  begin  to  diminish.  Triple¬ 
junction  cells  have  efficiencies  of  approximately  31%.  Quad- 
junction  cells  have  efficiencies  of  approximately  33%.  It  is 
expected  that  five-  and  six-junction  solar  cells  can  have 
efficiencies  of  approximately  35%  [36]. 

This  simulation  uses  Indium  Gallium  Nitride  as 
photovoltaic  material  in  a  simulation.  With  its  full- 
spectrum  characteristics,  the  efficiency  of  Indium  Gallium 
Nitride  has  the  potential  to  surpass  previous  results. 
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A. 


SINGLE-JUNCTION  SOLAR  CELL 


Figure  37  shows  a  graphical  representation  of  the 
single  junction  solar  cell.  The  thickness  of  the  emitter  is 
0.01  /urn  and  the  doping  is  lxl016/cm3.  The  thickness  of  the 
base  is  3  /urn  and  the  doping  is  lxl  016/ cm3. 

I - 

_ o 

Emitter  n+InGaN  0.01/xm  lel6cm 
Base  p+InGaN  3^m  le!6cm 


Figure  37.  Simple  single- j unction  InGaN  solar  cell. 

The  first  step  is  to  develop  single  junction  solar 
cells.  Optical  data  in  the  form  of  dielectric  constants  was 
found  in  [34]  for  In0.2oGao.8oN .  This  corresponds  to  a 
calculated  band  gap  of  2.66  eV.  Additional  optical  data  was 
found  in  [35]  for  Ino.57Gao.43N  (calculated  Eg=l  .  60  eV)  , 
Ino.68Gao.32N  (calculated  Eg=1.31  eV)  ,  and  Ino.7sGao.22N 
(calculated  Eg=l . 1  eV)  .  The  dielectric  constants  were  then 
converted  to  index  of  refraction  (n)  and  extinction 
coefficient  (k)  using  a  Matlab  script  (Appendix  B) .  The  band 
gap  and  the  optical  data  were  entered  into  the  input  deck. 
The  simulation  ran  in  Silvaco  Atlas.  From  the  log  file,  the 
current  and  voltage  data  was  extracted. 

In  the  single- j unction  case,  each  of  the  cells  had  AMO 
as  the  input  spectrum.  Figure  38  shows  the  IV  curves  of  four 
single-junction  solar  cells.  The  IV  curve  pattern 

corresponds  to  empirical  findings  of  other  multi j unction 
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solar  cells.  Typically,  the  cell  with  the  highest  voltage 
has  the  lowest  current.  The  cell  with  the  lowest  voltage  has 
the  highest  current. 


Solar  Cell  IV  Curve 


Table  7  shows  Isc,  Voc,  fill  factor,  and  efficiency  of 
the  single- j unction  solar  cells  presented  in  Fig.  36. 


Material 

Isc 

(mA/cm2) 

Voc  (V) 

Fill 

factor 

Efficiency 

Ino.20Gao.80N 

17 . 613 

2.23 

83.72% 

24 . 32% 

Ino.57Gao.43N 

18.395 

1  .  15 

78 . 80% 

12.35% 

Ino.68Gao.32N 

20.453 

0.91 

76.38% 

10.52% 

Ino.78Gao.22N 

20 . 848 

0.72 

73.73% 

08 .19% 

Table  7.  Efficiencies  of  four  single- j unction  InGaN  cells. 
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The  efficiencies  shown  in  Table  7  indicate  that  a  high- 
efficiency,  single- j  unction  cell  with  In0.2oGa0.8oN  has  good 
prospects.  A  24.32%  efficiency  compares  favorably  to  the 
single- j unction  Gallium  Arsenide  efficiency  of  25.1%.  The 
InGaN  cells  simulated  can  be  optimized  by  changing  the 
thickness  of  the  junction  and  changing  the  doping  levels. 
Adding  a  window,  a  back  surface  field  (BSF)  and  a  buffer  can 
also  help  in  improving  the  efficiency  results.  These 
improvements  can  be  the  subject  of  future  research. 


B.  DUAL- JUNCTION  SOLAR  CELL 

Figure  39  shows  a  graphical  representation  of  a  dual¬ 
junction  InGaN  solar  cell.  The  thicknesses  and  doping  levels 
remain  identical  to  the  single- j unction  case. 


From  the  discussion  in  the  single- j unction  solar  cell 
section,  it  should  be  clear  that  the  best  combination  of 
InGaN  band  gaps  is  2.66  eV  and  1.60  eV,  since  they  have  the 
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highest  efficiencies.  It  should  be  noted  that  for  this 
simulation,  the  AMO  spectrum  was  provided  for  the  top 
junction  In0.2oGao.8oN  (Eg=2.66  eV)  .  The  bottom  junction 
Ino.57Gao.43N  (Eg=l .  60  eV)  received  the  AMO  spectrum  minus  the 
spectrum  absorbed  by  the  top  cell. 

Figure  40  shows  the  IV  curves  of  the  dual- j unction 
solar  cell.  Note  that  compared  to  Figure  38,  the  IV  curve  of 
Eg=2.66  eV  stayed  the  same,  while  the  IV  curve  of  Eg=l . 60  eV 
had  a  current  drop.  This  current  decrease  was  expected 
because  the  input  spectrum  was  diminished  by  the  top  cell. 


Solar  Cell  IV  Curve 


Figure  40.  Dual- j unction  InGaN  solar  cell  IV  curve. 

Since  the  junctions  are  in  series,  the  overall  IV  curve 
is  limited  in  its  current  level  by  the  IV  curve  with  the 
lowest  current.  The  voltages  are  added  accordingly. 

Table  8  shows  the  efficiency  of  the  dual- j unction  solar 

cell . 
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Material 

Isc 

(mA/cm2) 

Voc  (V) 

Fill 

factor 

Efficiency 

Ino.2oGao.8oN 

In0.57Gao.43N 

Dual 

j  unction 

17 .11 

3.38 

83.25% 

35.58% 

Table  8.  Dual- j unction  InGaN  efficiency. 


There  is  a  significant  increase  in  efficiency  from 
single- j unction  (24.32%)  to  dual  junction  (35.58%)  InGaN 
solar  cells.  This  improvement  alone  should  be  enough  to 
continue  to  pursue  the  development  of  InGaN  solar  cells. 

Searching  for  further  improvements,  three-  and  four- 
junction  solar  cells  are  examined  next. 

C.  THREE- JUNCTION  SOLAR  CELL 

Figure  41  shows  a  graphical  representation  of  a  three- 
junction  InGaN  solar  cell.  Thicknesses  and  doping  levels 
stayed  the  same  as  the  single-  and  dual- j unction  cases.  From 
the  single- j unction  solar  cell  discussion,  the  three  most 
efficient  cells  were  selected  (Eg=2.66  eV,  Eg=l . 60  eV, 
Eg=1.31  eV)  .  The  top  cell  In0.2oGao.8oN  (Eg=2.66  eV)  received 
the  AMO  spectrum.  The  cell  Ino.57Gao.43N  (Eg=l .  60  eV)  received 
the  AMO  spectrum  minus  the  spectrum  absorbed  by  the  top 
cell.  The  cell  Ino.68Gao.32N  (Eg=1.31  eV)  received  the  AMO 
spectrum  minus  the  spectrum  absorbed  by  the  top  two  cells. 

Figure  42  shows  the  IV  curve  of  the  three- j unction 
solar  cell.  Note  that  the  current  level  of  the  top  two  cells 
is  identical  to  the  dual- j unction  case.  However,  the  current 
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level  of  the  bottom  cell  is  lower  than  the  single  junction 
case  for  the  same  concentration.  Similarly  to  the  dual¬ 
junction  case,  the  current  of  the  overall  IV  curve  is 
limited  by  the  individual  junction  with  the  lowest  IV  curve. 
The  voltages  are  added  accordingly,  since  the  three 
junctions  are  in  series. 
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Solar  Cell  IV  Curve 


Figure  42.  Three- j unction  InGaN  solar  cell  IV  curve. 


Table  9  shows  the  efficiency  of  the  three- j unction 
InGaN  solar  cell. 


Material 

Isc 

(mA/cm2) 

Voc  (V) 

Fill 

factor 

Efficiency 

Ino.20Gao.80N 

Ino.57Gao.43N 

Ino.  68Gao.  32N 

Three 

junction 

14 . 12 

4.27 

87.23% 

38 . 90% 

Table  9.  Three- j unction  InGaN  efficiency 


The  increase  in  efficiency  from  dual- j unction  to  three- 
junction  is  35.58%  to  38.90%.  The  modest  increase  in 
efficiency  can  be  attributed  to  the  drop  in  Isc  from  17.11 


56 


mA/cm2  to  14.12  mA/cm2 .  However,  the  increase  in  Voc  from 
3.38  V  to  4.27  ensured  that  the  efficiency  would  increase. 

The  last  case  to  be  examined  is  the  quad- j unction  InGaN 
solar  cell. 

D.  QUAD- JUNCTION  SOLAR  CELL 

Figure  43  shows  a  graphical  representation  of  a  quad- 
junction  InGaN  solar  cell.  Thicknesses  and  doping  levels 
remained  the  same  as  the  single-,  dual-,  and  three- j unction 
cases.  Spectrum  input  followed  the  pattern  of  the  dual-  and 
three- j unction  cases.  As  expected,  the  current  dropped  for 
all  cases  where  AMO  had  been  reduced  by  the  top  cells. 

Figure  44  shows  the  IV  curve  for  a  quad- j unction  InGaN 
solar  cell.  Note  that  the  IV-curves  for  the  top  three  cells 
are  identical  to  the  three- j unction  case.  The  bottom  cell 
has  a  lower  current  level  compared  to  the  single  junction 
case . 

Table  10  shows  the  efficiency  of  the  quad- j unction 
InGaN  solar  cell. 
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Solar  Cell  IV  Curve 


Figure  44.  Quad- j unction  InGaN  solar  cell  IV  curve 


Material 

Isc 

(mA/cm2) 

Voc  (V) 

Fill 

factor 

Efficiency 

Ino.20Gao.80N 

In0.57Gao.43N 

In0.68Gao.32N 

Ino.78Gao.22N 

Quad- 

junction 

12 . 88 

4 . 98 

87 .86% 

41 . 69% 

Table  10.  Quad 

junction  InGaN  efficiency 

For  completeness,  efficiency  calculations  are  provided 
below  for  the  quad- j unction  case.  These  calculations  were 
performed  with  the  help  of  a  Matlab  script. 
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Pmax= ImpVmp— (0.0 126 


-^)(4.4736  V)-0.0564 
cm' 


W 

cm2 


FF- 


P  max 
Isc  V  oc 


Imp  V mp 
Isc  V  oc 


0.0564— 

-  cm~ - -87.86% 

(0.0129— t)(4.9793V) 
cm 


P  max 


Imp  V  mp 

Pm 


0.0564 


W 


cm- 


0.1353 


W 


cm' 


-0.4169-41.69% 


In  order  to  get  a  perspective  on  how  this  result 
compares  to  actual  solar  cells  in  production.  Figure  45  is 
presented . 


Figure  45.  Spectrolab' s  solar  cell  efficiencies  (From  [34]). 
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The  label  nJ  stands  for  new  generation  four  to  six- 
junction  cells.  XTJ  is  a  designation  for  a  triple  junction 
cell.  UTJ  stands  for  Ultra  Triple  Junction.  ITJ  means 
Improved  Triple  Junction.  TJ  is  Triple  Junction.  DJ  is  Dual 
Junction.  SJ  is  single  junction. 

Figure  45  shows  the  efficiency  progression  of  a  leading 
solar  cell  manufacturer  (Spectrolab)  .  The  first  step 
consists  of  single  junction  silicon  solar  cells.  These  cells 
were  in  production  from  the  mid-1960s  until  the  early  1990s. 
Their  efficiency  was  close  to  15%.  The  next  step  consisted 
of  single  junction  gallium  arsenide  solar  cells.  These  cells 
had  efficiencies  between  15%  and  20%.  In  the  mid-1990s,  dual 
junction  solar  cells  appeared.  Efficiencies  were  slightly 
above  20%.  Triple  junction  cells  appeared  in  the  early  2000s 
with  efficiencies  nearing  30%.  The  projections  are  to 
produce  quad- j unction  solar  cells  within  the  next  decade 
with  efficiencies  just  under  35%. 

Further  comparisons  can  be  made  with  another 
simulation.  Modeling  in  [5]  used  physics  equations  to 
calculate  Isc  and  Voc.  The  input  spectrum  used  in  that 
simulation  was  AMI .  5  instead  of  AMO.  The  InGaN  band  gap 
formula  used  was  from  [37]  .  This  thesis  used  the  band  gap 
formula  from  [30]  .  Figure  4  6  shows  a  comparison  of  the  two 
formulas.  Both  formulas  seek  to  duplicate  findings  from 
actual  band  gap  measurements.  Therefore,  both  formulas  are 
approximations  only. 

The  simulation  from  [5]  obtained  the  results  presented 
in  Table  11.  The  highest  efficiency  obtained  in  that 
simulation  is  40.346%  with  a  six-junction  InGaN  solar  cell. 
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Indium  concentration  vs  InGaN  Band  Gap 


Figure  46.  Comparison  of  InGaN  band  gap  formulas 


Number  of 
junctions  in  the  cell 

/sc  (mA/cnr) 

Foe  (V) 

Output  peak  power 
(mW/cmj 

Efficiency  (%) 

2 

27 

1.22 

26.48 

27.485 

3 

18.2 

2.22 

32.42 

33.642 

4 

13.5 

3.28 

35.48 

36.825 

5 

11 

4.28 

37.68 

39.105 

6 

9.1 

5.34 

38.88 

40.346 

Table  1 1 . 

InGaN 

efficiency 

results  (From 

[5]  )  . 

The  differences  in  efficiency  results  can  be  partly 
attributed  to  the  difference  in  band  gap  formula  used.  For 
example,  In0.2oGao.8oN  has  a  band  gap  of  2.66  eV  with  the 
formula  used  in  this  thesis.  The  formula  from  [37]  yields  a 
band  gap  of  1.897  eV.  If  the  In0.2oGao.8oN  band  gap  of  1.897 
eV  is  substituted  in  the  single- j unction  solar  cell 
simulation  of  this  thesis,  the  efficiency  decreases 
significantly.  Figure  47  compares  the  IV  output  for 
Ino.20Gao.80N  at  Eg=1.897  eV  and  Eg=2.66  eV. 


62 


Solar  Cell  IV  Curve 


Figure  47.  IV  curve  for  Ino.20Gao.80N  using  different  band  gaps 


The  efficiency  of  the  Ino.20Gao.80N  (2.66  eV)  single¬ 
junction  solar  cell  simulated  in  this  thesis  was  24.32%. 
Changing  the  band  gap  to  1.897  eV  yields  an  efficiency  of 
14.94%.  This  is  a  significant  decrease.  However,  the  band 
gap  for  the  other  three  junctions  increases  with  the  formula 
used  in  [37] .  The  quad- j unction  simulation  was  run  with  the 
new  band  gaps.  Figure  4  8  shows  the  IV  curve  for  the  quad- 
junction  solar  cell.  The  efficiency  for  the  quad- j unction 
solar  cell  increases  to  43.62%.  It  is  difficult  to  state 
which  band  gap  model  is  more  correct. 
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Solar  Cell  IV  Curve 


Figure  48.  Quad- j unction  InGaN  solar  cell  IV  curve  using 
calculated  band  gaps  from  [37]  formula. 


Material 

Isc 

(mA/cm2) 

Voc  (V) 

Fill 

factor 

Efficiency 

Ino.2oGao.8oN 

In0.57Gao.43N 

Ino.68Gao.32N 

In0.78Gao.22N 

Quad- 

junction 

12 . 9 

5.194 

88.06% 

43.62 

Table  12.  Quad  junction  InGaN  efficiency  using  band  gaps 

from  [37]  calculations. 
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A  common  finding  in  [5]  and  this  thesis  is  that  both 
simulations  show  that  Indium  Gallium  Nitride  multi j unction 
solar  cells  can  provide  a  significant  improvement  in  solar 
cell  efficiency.  However,  this  is  the  first  time  that  a 
Technology  Computer  Aided  Design  (TCAD) ,  such  as  Silvaco 
Atlas,  has  been  used  to  simulate  InGaN  solar  cells. 

The  results  of  this  thesis  simulation  demonstrate  that 
Indium  Gallium  Nitride  is  potentially  an  excellent 
semiconductor  photovoltaic  material.  The  material  science 
research  in  the  future  can  confirm  these  outcomes  in  the 
future . 

Production  of  single- j unction  solar  cells  appears  to  be 
the  first  step  to  be  taken.  Once  the  photovoltaic  properties 
are  demonstrated  with  actual  Indium  Gallium  Nitride,  the 
complexity  of  dual- j unction  solar  cells  can  be  addressed. 
Tunnel  junctions  made  up  of  Indium  Gallium  Nitride  may  need 
to  be  explored  as  well. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  RESULTS  AND  CONCLUSIONS 

Indium  Gallium  Nitride  is  a  semiconductor  material  with 
potential  to  be  used  in  photovoltaic  devices.  A  new 
simulation  was  performed  using  Silvaco  Atlas.  The  results  of 
the  quad- j unction  Indium  Gallium  Nitride  solar  cell  indicate 
that  a  new  high-efficiency  material  should  be  produced. 
Current  state-of-the-art  multi j unction  solar  cells  have 
efficiencies  in  the  30-33%  range.  The  41%  efficiency 
predicted  in  this  thesis  is  the  highest  of  all  simulations 
performed  at  the  Naval  Postgraduate  School  using  Silvaco 
Atlas . 

B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

There  are  multiple  areas  that  can  be  explored  in  future 
research.  This  thesis  focused  on  the  development  of  a  solar 
cell  model  that  emphasizes  the  use  of  optical  constants 
(refraction  and  extinction  coefficients  n  and  k) .  The  model 
used  default  settings  for  other  parameters,  such  as 
permittivity,  affinity,  radiative  recombination  rate, 
electron  and  hole  lifetimes,  electron  and  hole  density  of 
states,  and  lattice  constants.  One  area  of  future  research 
is  to  obtain  measured  data  for  the  above  parameters.  This 
subject  is  critical  in  improving  the  model.  The  other  option 
is  to  interpolate  from  known  values.  Additionally,  the 
inclusion  of  tunnel  junctions  in  the  simulation  is  of  great 
value . 
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Other  areas  of  research  include  finding  different 
semiconductor  materials  that  offer  potential  for  high- 
efficiency  solar  cells.  Lawrence  Berkeley  National 
Laboratory  is  also  working  with  multiband  materials,  such  as 
Zinc  Manganese  Tellerium  Oxide  (ZnMnTeO) .  This  type  of 
material  offers  a  single  junction  solar  cell,  but  there  are 
multiple  band  gaps  within  the  single  junction. 

One  more  recommended  path  is  to  optimize  the  physical 
parameters  of  the  solar  cell.  The  thickness  of  each  of  the 
layers  can  be  optimized  to  produce  a  better  efficiency.  The 
band  gap  distribution  for  a  quad- j unction  solar  cell  can  be 
improved  to  obtain  higher  efficiency.  It  should  be  noted 
that  obtaining  the  correct  optical  data  for  these  band  gaps 
is  essential. 
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APPENDIX  A:  SILVACO  ATLAS  INPUT  DECK 


The  code  for  the  quad- j  unction  InGaN  solar  cell  was 
broken  down  into  four  single- j unction  solar  cells.  This  was 
done  because  the  material  properties  of  InGaN  had  to  be 
changed  for  each  of  the  junctions. 

Original  code  for  the  Silvaco  Atlas  input  deck  was 
obtained  from  Michalopolous  [1],  Bates,  [2],  Green  [3],  and 
Canfield  [4] .  Modifications  to  the  code  were  made  to  support 
this  thesis. 

A.  TOP  JUNCTION:  IN0.2oGA0.8oN,  EG=2 . 66  EV 

go  atlas 


set  cellWidth=5 . 000000e+002 

set  capWidthpercent=8 . 000000e+000 

set  divs=l . 000000e+001 

set  contThick=l . 000000e-001 

set  capThick=3 . 000000e-001 

set  capDop=l . 000000e+020 

set  windowThick=0 . 01 

set  winDop=2 . 15el7 

set  emitterThick=0 . 01 

#changed  emitDop  from  lel6  to  le20 

set  emitDop=lel6 

set  baseThick=3 . 1 94 67 

#changed  basDop  from  lel6  to  le20 

set  baseDop=lel6 

set  bsfThick=0 . 03533 

set  bsf Dop=2 . 15el 9 

set  cellWidthDiv=$cellWidth/ $divs 
set  width3d=100e6/$cellWidth 

set  capWidth=0 . 01*$capWidthpercent*$cellWidth/2 
set  capWidthDiv=$capWidth/ ($divs/2) 
set  cellWidthHalf =$cellWidth/2 

set  bsfLo=0 

set  bsfHi=$bsfLo-$bsfThick 
set  bsf Div=$bsf Thick/ $divs 

set  baseLo=$bsfHi 
set  baseHi=$baseLo-$baseThick 
set  baseMid=$baseLo-$baseThick/2 
set  baseDiv=$baseThick/ $divs 

set  emitterLo=$baseHi 
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set  emi t t er Hi =$emitterLo-$emit ter Thick 
set  emitterDiv=$emitterThick/ $divs 

set  windowLo=$emitterHi 

set  windowHi=$windowLo-$windowThick 

set  windowDiv=$windowThick/ $divs 

set  capLo=$windowHi 

set  capHi=$capLo-$capThick 

#set  capDiv=$capThick/$divs 

set  contLo=$capHi 

set  contHi=$contLo-$contThick 

set  contDiv=$contThick/ $divs 

set  lightY=$emitterHi-5 

mesh  width=$width3d 
##  X-Mesh 

x.mesh  loc=-$cellWidthHalf  spac=$cellWidthDiv 
x.mesh  loc=-$capWidth  spac=$capWidthDiv 
x.mesh  loc=$capWidth  spac=$capWidthDiv 

x. mesh  loc=$cellWidthHalf  spac=$cellWidthDiv 

##  Y-Mesh 

#  Top  contact 

y. mesh  loc=$contHi  spac=0 
y.mesh  loc=$contLo  spac=0 

#  Cap 

#  Window 

y.mesh  loc=$windowHi  spac=$windowDiv 
y.mesh  loc=$windowLo  spac=$windowDiv 

#  Emitter 

y.mesh  loc=$emitterLo  spac=$emitterDiv 

#  Base 

y.mesh  loc=$baseMid  spac=$baseDiv 

#  BSF 

y.mesh  loc=$bsfHi  spac=$bsfDiv 
y.mesh  loc=$bsfLo  spac=$bsfDiv 


##  Regions 

#  Cap 

#region  num=8  material=Vacuum  x .min=-$capWidth  x . max=$capWidth  y.min=$contHi 
y ,max=$contLo 

region  num=l  material=InGaN  x . min=-$capWidth  x . max=$capWidth  y.min=$capHi 
y.max=$capLo  x.comp=0.20 

region  num=2  material=Vacuum  x .min=-$cellWidthHalf  x ,max=-$capWidth 
y ,min=$contHi  y.max=$capLo 

region  num=3  material=Vacuum  x ,min=$capWidth  x.max=$cellWidthHalf  y.min=$contHi 
y ,max=$capLo 

#  Window  [for  Ge  cell,  use  AlGaAs  with  x.comp=0.2] 

region  num=4  material=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowLo  x.comp=0.20 


#region  num=4  material=AlGaAs  x.comp=0.2  x.min=-$cellWidthHalf 

x. max=$cellWidthHalf  y.min=$windowHi  y.max=$windowLo 

#  Emitter 

region  num=5  mater ial=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y . min=$emitterHi  y.max=$emitterLo  x.comp=0.20 

#  Base 
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region  num=6  mater ial=InGaN  x .min=-$cellWidthHalf  x .max=$cellWidthHalf 
y ,min=$baseHi  y ,max=$baseLo  x.comp=0.20 

♦  BSF 

region  num=7  mater ial=InGaN  x .min=-$cellWidthHalf  x .max=$cellWidthHalf 
y.min=$bsfHi  y.max=$bsfLo  x.comp=0.20 

##  Electrodes  [for  InGaP  cell,  add  cathode  (gold)  and  remove  cathode 
(conductor) ] 

electrode  name=cathode  material=Gold  x  .min=-$capWidth  x . max=$capWidth 
y ,min=$contHi  y ,max=$contLo 

#electrode  name=cathode  x .min=-$cellWidthHalf  x .max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowHi 

electrode  name=anode  x.min=-$cellWidthHalf  x .max=$cellWidthHalf  y.min=$bsfLo 
y ,max=$bsfLo 

##  Doping  [for  InGaP  cell,  uncomment  cap  doping] 

♦  Cap 

doping  uniform  region=l  n.type  conc=$capDop 

♦  Window 

doping  uniform  region=4  n.type  conc=$winDop 

♦  Emitter 

doping  uniform  region=5  n.type  conc=$emitDop 

♦  Base 

doping  uniform  region=6  p.type  conc=$baseDop 

♦  BSF 

doping  uniform  region=7  p.type  conc=$bsfDop 
##  Material  properties 

♦  Opaque  contact  [comment  out  for  InGaP  cell] 

♦material  region=8  real . index=l . 2  imag . index=l . 8 

♦  Vacuum  (for  zero  reflection)  [change  to  match  window  material  (InGaP  use 
Vacuum_AlInP) ] 

♦  [for  InGaP  cell,  comment  out  region  1] 

♦material  region=l  index. file=Vacuum_InGaP. opt 

material  region=2  index . f ile=VacuumIn20Ga80N . opt 
material  region=3  index . f ile=VacuumIn20Ga80N . opt 

♦inGaN 

material  material=InGaN  EG300=2.6612  index . file=In2 0Ga8 ON . opt 

♦  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 


♦♦  Models  [InGaP  cell,  1;  GaAs  cell,  5&6;  InGaNAs  cell,  7] 
models  region=l  CONMOB 

MODELS  CHUANG  CONMOB  FLDMOB  SRH  OPTR  PRINT 

♦  ♦  Light  beams  [GaAs  bl ,  0 . 55-0 . 75 , 2 00  b2 ,  0 . 75-0 . 8 8 ,  65 ]  0.12-2.7,50  [630,825] 
beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . file=AM0nrel . spec  \ 

wavel . start=0 . 12  wavel . end=2 . 4  wavel.num=50 

struct  outf ile=SingleCell_webf . str 
♦tonyplot  SingleCell_webf . str 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=0 . 9 

♦♦  Getting  Isc  for  I-V  curve  points 
method  newton  maxtraps=10  itlimit=100 
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solve  bl=0 . 95 

extract  name="isc"  max ( i . "cathode" ) 

set  isc=$isc*$width3d 

set  il=$isc/10 

set  i2=$il+$isc/10 

set  i3=$i2+$isc/10 

set  i4=$i3+$isc/10 

set  i5=$i4+$isc/10 

set  i6=$i5+$isc/20 

set  i7=$i6t$isc/20 

set  i8=$i7+$isc/20 

set  i9=$i8+$isc/20 

set  il0=$i9+$isc/20 

set  ill=$il0+$isc/40 

set  il2=$ill+$isc/40 

set  il3=$il2+$isc/40 

set  il4=$il3+$isc/40 

set  il5=$il4+$isc/40 

set  il 6=$il 5+$isc/8 0 

set  il7=$il 6+$isc/8 0 

set  il8=$il 7+$isc/8 0 

set  il9=$il8+$isc/80 

set  i20=$il 9+$isc/8 0 

set  i21=$i2 0+$isc/8 0 

set  i22=$i2 l+$isc/8 0 

set  i23=$i22+$isc/8 0 

set  i24=$i23+$isc/8 0 

set  i25=$i24+$isc/80-0. 00001 

## 

log  outfile=In20Ga80N . log 

method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 

contact  name=anode  current 

method  newton  maxtraps=10  itlimit=100 

##  Pmax  points  [InGaP  18-25;  GaAs  15-25;  InGaNAs  13-25;  Ge  11-25] 


solve 

ianode= 

-$i25 

bl=0 

.  95 
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-$i24 

bl=0 
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-$i2 1 
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ianode= 
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95 
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95 
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95 
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solve  ianode=-$i2  bl=0.95 
solve  ianode=-$il  bl=0.95 

solve  ianode=0  bl=0.95 
log  off 

extract  name="iv"  curve  (v . "anode" ,  i. "cathode")  outf ile=" IVcurveIn2 0Ga8 ON . dat " 
tonyplot  IVcurveIn20Ga80N.dat 

log  outf ile=doneIn2 0Ga8 ON . log 
log  off 

B.  SECOND  JUNCTION:  IN0.57GA0.43N,  EG=1 . 6  EV 

go  atlas 

set  cellWidth=5 . 000000e+002 

set  capWidthpercent=8 . OOOOOOe+OOO 

set  divs=l . OOOOOOe+OOl 

set  contThick=l . 000000e-001 

set  capThick=3 . 000000e-001 

set  capDop=l . 000000e+020 

set  windowThick=0 . 01 

set  winDop=2 . 15el7 

set  emitterThick=0 . 01 

#changed  emitDop  from  lel6  to  le20 

set  emitDop=lel6 

set  baseThick=3 . 1 94 67 

#changed  basDop  from  lel6  to  le20 

set  baseDop=lel6 

set  bsfThick=0 . 03533 

set  bsf Dop=2 . 15el 9 

set  cellWidthDiv=$cellWidth/ $divs 
set  width3d=100e6/$cellWidth 

set  capWidth=0 . 01*$capWidthpercent*$cellWidth/2 
set  capWidthDiv=$capWidth/ ($divs/2) 
set  cellWidthHalf =$cellWidth/2 

set  bsfLo=0 

set  bsfHi=$bsfLo-$bsfThick 
set  bsf Div=$bsf Thick/ $divs 

set  baseLo=$bsfHi 
set  baseHi=$baseLo-$baseThick 
set  baseMid=$baseLo-$baseThick/2 
set  baseDiv=$baseThick/ $divs 

set  emitterLo=$baseHi 

set  emi t t er Hi =$emitterLo-$emit ter Thick 
set  emitterDiv=$emitterThick/ $divs 

set  windowLo=$emitterHi 

set  windowHi=$windowLo-$windowThick 

set  windowDiv=$windowThick/ $divs 

set  capLo=$windowHi 

set  capHi=$capLo-$capThick 

#set  capDiv=$capThick/$divs 

set  contLo=$capHi 

set  contHi=$contLo-$contThick 
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set  contDiv=$contThick/ $divs 


set  lightY=$emitterHi-5 

mesh  width=$width3d 
##  X-Mesh 

x.mesh  loc=-$cellWidthHalf  spac=$cellWidthDiv 
x.mesh  loc=-$capWidth  spac=$capWidthDiv 
x.mesh  loc=$capWidth  spac=$capWidthDiv 

x. mesh  loc=$cellWidthHalf  spac=$cellWidthDiv 

##  Y-Mesh 

#  Top  contact 

y. mesh  loc=$contHi  spac=0 
y.mesh  loc=$contLo  spac=0 

#  Cap 

#  Window 

y.mesh  loc=$windowHi  spac=$windowDiv 
y.mesh  loc=$windowLo  spac=$windowDiv 

#  Emitter 

y.mesh  loc=$emitterLo  spac=$emitterDiv 

#  Base 

y.mesh  loc=$baseMid  spac=$baseDiv 

#  BSF 

y.mesh  loc=$bsfHi  spac=$bsfDiv 
y.mesh  loc=$bsfLo  spac=$bsfDiv 


##  Regions 

#  Cap 

#region  num=8  material=Vacuum  x ,min=-$capWidth  x . max=$capWidth  y.min=$contHi 
y ,max=$contLo 

region  num=l  material=InGaN  x . min=-$capWidth  x . max=$capWidth  y.min=$capHi 
y.max=$capLo  x.comp=0.57 

region  num=2  material=Vacuum  x ,min=-$cellWidthHalf  x ,max=-$capWidth 
y ,min=$contHi  y.max=$capLo 

region  num=3  material=Vacuum  x ,min=$capWidth  x.max=$cellWidthHalf  y.min=$contHi 
y ,max=$capLo 

#  Window  [for  Ge  cell,  use  AlGaAs  with  x.comp=0.2] 

region  num=4  material=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowLo  x.comp=0.57 


#region  num=4  material=AlGaAs  x.comp=0.2 

x. max=$cellWidthHalf  y.min=$windowHi  y.max=$windowLo 

#  Emitter 

region  num=5  mater ial=InGaN  x ,min=-$cellWidthHalf 

y , min=$emitterHi  y.max=$emitterLo  x.comp=0.57 

#  Base 

region  num=6  mater ial=InGaN  x ,min=-$cellWidthHalf 

y ,min=$baseHi  y ,max=$baseLo  x.comp=0.57 

#  BSF 

region  num=7  mater ial=InGaN  x ,min=-$cellWidthHalf 

y.min=$bsfHi  y.max=$bsfLo  x.comp=0.57 


x . min=-$cellWidthHalf 


x ,max=$cellWidthHalf 


x ,max=$cellWidthHalf 


x ,max=$cellWidthHalf 


##  Electrodes  [for  InGaP  cell,  add  cathode  (gold)  and  remove  cathode 
(conductor) ] 

electrode  name=cathode  material=Gold  x ,min=-$capWidth  x . max=$capWidth 
y ,min=$contHi  y ,max=$contLo 

#electrode  name=cathode  x ,min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowHi 

electrode  name=anode  x.min=-$cellWidthHalf  x ,max=$cellWidthHalf  y.min=$bsfLo 
y ,max=$bsfLo 
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##  Doping  [for 

#  Cap 

doping  uniform 

#  Window 
doping  uniform 

#  Emitter 
doping  uniform 

#  Base 

doping  uniform 

#  BSF 

doping  uniform 


InGaP  cell,  uncomment  cap  doping] 
region=l  n.type  conc=$capDop 
region=4  n.type  conc=$winDop 
region=5  n.type  conc=$emitDop 
region=6  p.type  conc=$baseDop 
region=7  p.type  conc=$bsfDop 


##  Material  properties 

#  Opaque  contact  [comment  out  for  InGaP  cell] 

♦material  region=8  real . index=l . 2  imag . index=l . 8 

#  Vacuum  (for  zero  reflection)  [change  to  match  window  material  (InGaP  use 
Vacuum_AlInP) ] 

#  [for  InGaP  cell,  comment  out  region  1] 
material  region=2  index . f ile=VacuumIn57Ga43N . opt 

material  region=3  index . f ile=VacuumIn57Ga43N . opt 

#  GaAs 


#InGaN 

material  material=InGaN  EG300=1.559  index . file=In57Ga43N . opt 
#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 


##  Models  [InGaP  cell,  1;  GaAs  cell,  5&6;  InGaNAs  cell,  7] 
models  region=l  CONMOB 

MODELS  CHUANG  CONMOB  FLDMOB  SRH  OPTR  PRINT 

##  Light  beams  [GaAs  bl , 0 . 55-0 . 75 , 2 00  b2 , 0 . 75-0 . 8 8 , 65 ]  0.12-2.7,50  [630,825] 
beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . f ile=Post Junctionl . spec . txt  \ 

wavel . start=0 . 12  wavel . end=2 . 4  wavel.num=50 

struct  outf ile=SingleCell_webf . str 
♦tonyplot  SingleCell_webf . str 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=0 . 9 

##  Getting  Isc  for  I-V  curve  points 
method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 

extract  name="isc"  max ( i . "cathode" ) 

set  isc=$isc*$width3d 

set  il=$isc/10 

set  i2=$il+$isc/10 

set  i3=$i2+$isc/10 

set  i4=$i3t$isc/10 

set  i5=$i4+$isc/10 

set  i6=$i5t$isc/20 

set  i7=$i6t$isc/20 

set  i8=$i7+$isc/20 

set  i9=$i8t$isc/20 

set  il0=$i9t$isc/20 
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set  ill=$il0+$isc/40 

set  il2=$ill+$isc/40 

set  il3=$il2+$isc/40 

set  il4=$il3+$isc/40 

set  il5=$il4+$isc/40 

set  il 6=$il 5+$isc/8 0 

set  il7=$il6+$isc/80 

set  il8=$il 7+$isc/8 0 

set  il9=$il8+$isc/80 

set  i20=$il 9+$isc/8 0 

set  i21=$i2 0+$isc/8 0 

set  i22=$i2 l+$isc/8 0 

set  i23=$i22+$isc/8 0 

set  i24=$i23+$isc/8 0 

set  i25=$i24+$isc/80-0. 00001 

## 

log  outfile=In57Ga43N . log 

method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 

contact  name=anode  current 

method  newton  maxtraps=10  itlimit=100 

##  Pmax  points  [InGaP  18-25;  GaAs  15-25;  InGaNAs  13-25;  Ge  11-25] 


solve 

ianode= 

-$i2  5 

bl=0 

.  95 

solve 

ianode= 

-$i24 

bl=0 

.  95 

solve 

ianode= 

-$i23 

bl=0 

.  95 

solve 

ianode= 

-$i22 

bl=0 

.  95 

solve 

ianode= 

-$i2 1 

bl=0 

.  95 

solve 

ianode= 

-$i20 

bl=0 

.  95 

solve 

ianode= 

-$il  9 

bl=0 

.  95 

solve 

ianode= 

-$il8 

bl=0 

.  95 

solve 

ianode= 

-$il7 

bl=0 

.  95 

solve 

ianode= 

-$il  6 

bl=0 

.  95 

solve 

ianode= 

-$il5 

bl=0 

.  95 

solve 

ianode= 

-$il4 

bl=0 

.  95 

solve 

ianode= 

-$il3 

bl=0 

.  95 

solve 

ianode= 

-$il2 

bl=0 

.  95 

solve 

ianode= 

-$ill 

bl=0 

.  95 

solve 

ianode= 

-$il0 

bl=0 

.  95 

solve 

ianode= 

-$19 

bl=0 . 

95 

solve 

ianode= 

-$i8 

bl  =  0 . 

95 

solve 

ianode= 

-$17 

bl  =  0 . 

95 

solve 

ianode= 

-$i6 

bl  =  0 . 

95 

solve 

ianode= 

-$15 

bl  =  0  . 

95 

solve 

ianode= 

-$14 

bl  =  0 . 

95 

solve 

ianode= 

-$13 

bl  =  0 . 

95 

solve 

ianode= 

-$12 

bl  =  0 . 

95 

solve 

ianode= 

-$il 

bl  =  0 . 

95 

## 

solve  ianode=0  bl=0.95 
log  off 

extract  name="iv"  curve  (v . "anode" ,  i. "cathode")  outf ile=" IVcurveIn57Ga4 3N . dat " 
tonyplot  IVcurveIn57Ga43N.dat 


log  outf He=doneIn57Ga4 3N . log 
log  off 
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C.  THIRD  JUNCTION:  IN0.6sGA0.32N,  EG=1.31  EV 

go  atlas 


set  cellWidth=5 . 000000e+002 

set  capWidthpercent=8 . 000000e+000 

set  divs=l . 000000e+001 

set  contThick=l . 000000e-001 

set  capThick=3 . 000000e-001 

set  capDop=l . 000000e+020 

set  windowThick=0 . 01 

set  winDop=2 . 15el7 

set  emitterThick=0 . 01 

#changed  emitDop  from  lel6  to  le20 

set  emitDop=lel6 

set  baseThick=3 . 1 94 67 

#changed  basDop  from  lel6  to  le20 

set  baseDop=lel6 

set  bsfThick=0 . 03533 

set  bsf Dop=2 . 15el 9 

set  cellWidthDiv=$cellWidth/ $divs 
set  width3d=100e6/$cellWidth 

set  capWidth=0 . 01*$capWidthpercent*$cellWidth/2 
set  capWidthDiv=$capWidth/ ($divs/2) 
set  cellWidthHalf =$cellWidth/2 

set  bsfLo=0 

set  bsfHi=$bsfLo-$bsfThick 
set  bsf Div=$bsf Thick/ $divs 

set  baseLo=$bsfHi 
set  baseHi=$baseLo-$baseThick 
set  baseMid=$baseLo-$baseThick/2 
set  baseDiv=$baseThick/ $divs 

set  emitterLo=$baseHi 

set  emi t t er Hi =$emitterLo-$emit ter Thick 
set  emitterDiv=$emitterThick/ $divs 

set  windowLo=$emitterHi 

set  windowHi=$windowLo-$windowThick 

set  windowDiv=$windowThick/ $divs 

set  capLo=$windowHi 

set  capHi=$capLo-$capThick 

#set  capDiv=$capThick/$divs 

set  contLo=$capHi 

set  contHi=$contLo-$contThick 

set  contDiv=$contThick/ $divs 

set  lightY=$emitterHi-5 

mesh  width=$width3d 
##  X-Mesh 

x.mesh  loc=-$cellWidthHalf  spac=$cellWidthDiv 
x.mesh  loc=-$capWidth  spac=$capWidthDiv 
x.mesh  loc=$capWidth  spac=$capWidthDiv 
x.mesh  loc=$cellWidthHalf  spac=$cellWidthDiv 
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##  Y-Mesh 

#  Top  contact 

y.mesh  loc=$contHi  spac=0 
y.mesh  loc=$contLo  spac=0 

#  Cap 

#  Window 

y.mesh  loc=$windowHi  spac=$windowDiv 
y.mesh  loc=$windowLo  spac=$windowDiv 

#  Emitter 

y.mesh  loc=$emitterLo  spac=$emitterDiv 

#  Base 

y.mesh  loc=$baseMid  spac=$baseDiv 

#  BSF 

y.mesh  loc=$bsfHi  spac=$bsfDiv 
y.mesh  loc=$bsfLo  spac=$bsfDiv 


##  Regions 

#  Cap 

#region  num=8  mater ial=Vacuum  x ,min=-$capWidth  x . max=$capWidth  y.min=$contHi 
y ,max=$contLo 

region  num=l  material=InGaN  x . min=-$capWidth  x . max=$capWidth  y.min=$capHi 
y.max=$capLo  x.comp=0.68 

region  num=2  material=Vacuum  x ,min=-$cellWidthHalf  x ,max=-$capWidth 
y ,min=$contHi  y.max=$capLo 

region  num=3  material=Vacuum  x ,min=$capWidth  x.max=$cellWidthHalf  y.min=$contHi 
y ,max=$capLo 

#  Window  [for  Ge  cell,  use  AlGaAs  with  x.comp=0.2] 

region  num=4  mater ial=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowLo  x.comp=0.68 

#region  num=4  material=AlGaAs  x.comp=0.2  x.min=-$cellWidthHalf 

x. max=$cellWidthHalf  y.min=$windowHi  y.max=$windowLo 

#  Emitter 

region  num=5  mater ial=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y . min=$emitterHi  y.max=$emitterLo  x.comp=0.68 

#  Base 

region  num=6  mater ial=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$baseHi  y ,max=$baseLo  x.comp=0.68 

#  BSF 

region  num=7  mater ial=InGaN  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y.min=$bsfHi  y.max=$bsfLo  x.comp=0.68 

##  Electrodes  [for  InGaP  cell,  add  cathode  (gold)  and  remove  cathode 

(conductor) ] 

electrode  name=cathode  material=Gold  x ,min=-$capWidth  x . max=$capWidth 
y ,min=$contHi  y ,max=$contLo 

#electrode  name=cathode  x .min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowHi 

electrode  name=anode  x.min=-$cellWidthHalf  x ,max=$cellWidthHalf  y.min=$bsfLo 
y ,max=$bsfLo 

##  Doping  [for  InGaP  cell,  uncomment  cap  doping] 

#  Cap 

doping  uniform  region=l  n.type  conc=$capDop 

#  Window 

doping  uniform  region=4  n.type  conc=$winDop 

#  Emitter 

doping  uniform  region=5  n.type  conc=$emitDop 

#  Base 

doping  uniform  region=6  p.type  conc=$baseDop 
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#  BSF 

doping  uniform  region=7  p.type  conc=$bsfDop 
##  Material  properties 

#  Opaque  contact  [comment  out  for  InGaP  cell] 

♦material  region=8  real . index=l . 2  imag . index=l . 8 

#  Vacuum  (for  zero  reflection) 

material  region=2  index . f ile=VacuumIn68Ga32N . opt 
material  region=3  index . f ile=VacuumIn68Ga32N . opt 

#InGaN 

material  material=InGaN  EG300=1.3068  index . file=In68Ga32N . opt 

#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 


##  Models  [InGaP  cell,  1;  GaAs  cell,  5&6;  InGaNAs  cell,  7] 

models  region=l  CONMOB 

MODELS  CONMOB  FLDMOB  SRH  OPTR  PRINT 

##  Light  beams  [GaAs  bl , 0 . 55-0 . 75 , 2 00  b2 , 0 . 75-0 . 8 8 , 65 ]  0.12-2.7,50  [630,825] 
beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . f ile=Post Junction2 . spec . txt  \ 

wavel . start=0 . 12  wavel . end=2 . 4  wavel.num=50 

struct  outf ile=SingleCell_webf . str 
♦tonyplot  SingleCell_webf . str 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=0 . 9 

##  Getting  Isc  for  I-V  curve  points 
method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 

extract  name="isc"  max ( i . "cathode" ) 

set  isc=$isc*$width3d 

set  il=$isc/10 

set  i2=$il+$isc/10 

set  i3=$i2+$isc/10 

set  i4=$i3+$isc/10 

set  i5=$i4+$isc/10 

set  i6=$i5+$isc/20 

set  i7=$i6+$isc/20 

set  i8=$i7+$isc/20 

set  i9=$i8+$isc/20 

set  il0=$i9+$isc/20 

set  ill=$il0+$isc/40 

set  il2=$ill+$isc/40 

set  il3=$il2+$isc/40 

set  il4=$il3+$isc/40 

set  il5=$il4+$isc/40 

set  il 6=$il 5+$isc/8 0 

set  il7=$il6+$isc/80 

set  il8=$il 7+$isc/8 0 

set  il9=$il8+$isc/80 

set  i20=$il 9+$isc/8 0 

set  i21=$i2 0+$isc/8 0 

set  i22=$i2 l+$isc/8 0 

set  i23=$i22+$isc/8 0 

set  i24=$i23+$isc/8 0 
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set  i25=$i24+$isc/80-0. 00001 

## 


log  outfile=In68Ga32N . log 

method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 

contact  name=anode  current 

method  newton  maxtraps=10  itlimit=100 

##  Pmax  points  [InGaP  18-25;  GaAs  15-25;  InGaNAs  13-25;  Ge  11-25] 


solve 

ianode= 

-$i25 

bl=0 

.  95 

solve 

ianode= 

-$i24 

bl=0 

.  95 

solve 

ianode= 

-$i23 

bl=0 

.  95 

solve 

ianode= 

-$i22 

bl=0 

.  95 

solve 

ianode= 

-$i2 1 

bl=0 

.  95 

solve 

ianode= 

-$i20 

bl=0 

.  95 

solve 

ianode= 

-$il  9 

bl=0 

.  95 

solve 

ianode= 

-$il8 

bl=0 

.  95 

solve 

ianode= 

-$il7 

bl=0 

.  95 

solve 

ianode= 

-$il  6 

bl=0 

.  95 

solve 

ianode= 

-$il5 

bl=0 

.  95 

solve 

ianode= 

-$il4 

bl=0 

.  95 

solve 

ianode= 

-$il3 

bl=0 

.  95 

solve 

ianode= 

-$il2 

bl=0 

.  95 

solve 

ianode= 

-$ill 

bl=0 

.  95 

solve 

ianode= 

-$il0 

bl=0 

.  95 

solve 

ianode= 

-$19 

bl=0 . 

95 

solve 

ianode= 

-$i8 

bl  =  0 . 

95 

solve 

ianode= 

-$17 

bl  =  0 . 

95 

solve 

ianode= 

-$i6 

bl  =  0 . 

95 

solve 

ianode= 

-$15 

bl  =  0 . 

95 

solve 

ianode= 

-$14 

bl  =  0 . 

95 

solve 

ianode= 

-$13 

bl  =  0 . 

95 

solve 

ianode= 

-$12 

bl  =  0 . 

95 

solve 

ianode= 

-$il 

bl  =  0 . 

95 

## 

solve  ianode=0  bl=0.95 
log  off 

extract  name="iv"  curve (v . "anode" ,  i. "cathode")  outf ile=" IVcurveIn68Ga32N . dat " 
tonyplot  IVcurveIn68Ga32N.dat 

log  outfile=doneIn68Ga32N . log 
log  off 


D.  BOTTOM  JUNCTION:  IN0.7sGA0.22N,  EG=1.11  EV 

go  atlas 


set  cellWidth=5 . 000000e+002 

set  capWidthpercent=8 . 000000e+000 

set  divs=l . 000000e+001 

set  contThick=l . 000000e-001 

set  capThick=3 . 000000e-001 

set  capDop=l . 000000e+020 
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set  windowThick=0 . 01 

set  winDop=2 . 15el7 

set  emitterThick=0 . 01 

#changed  emitDop  from  lel6  to  le20 

set  emitDop=lel6 

set  baseThick=3 . 1 94 67 

#changed  basDop  from  lel6  to  le20 

set  baseDop=lel6 

set  bsfThick=0 . 03533 

set  bsf Dop=2 . 15el 9 

set  cellWidthDiv=$cellWidth/$divs 
set  width3d=100e6/$cellWidth 

set  capWidth=0 . 01*$capWidthpercent*$cellWidth/2 
set  capWidthDiv=$capWidth/ ($divs/2) 
set  cellWidthHalf =$cellWidth/2 

set  bsfLo=0 

set  bsfHi=$bsfLo-$bsfThick 
set  bsf Div=$bsf Thick/ $divs 

set  baseLo=$bsfHi 
set  baseHi=$baseLo-$baseThick 
set  baseMid=$baseLo-$baseThick/2 
set  baseDiv=$baseThick/$divs 

set  emitterLo=$baseHi 

set  emi t t er Hi =$emitterLo-$emit ter Thick 
set  emitterDiv=$emitterThick/ $divs 

set  windowLo=$emitterHi 

set  windowHi=$windowLo-$windowThick 

set  windowDiv=$windowThick/ $divs 

set  capLo=$windowHi 

set  capHi=$capLo-$capThick 

#set  capDiv=$capThick/$divs 

set  contLo=$capHi 

set  contHi=$contLo-$contThick 

set  contDiv=$contThick/$divs 

set  lightY=$emitterHi-5 

mesh  width=$width3d 
##  X-Mesh 

x.mesh  loc=-$cellWidthHalf  spac=$cellWidthDiv 
x.mesh  loc=-$capWidth  spac=$capWidthDiv 
x.mesh  loc=$capWidth  spac=$capWidthDiv 

x. mesh  loc=$cellWidthHalf  spac=$cellWidthDiv 

##  Y-Mesh 

#  Top  contact 

y. mesh  loc=$contHi  spac=0 
y.mesh  loc=$contLo  spac=0 

#  Cap 

#  Window 

y.mesh  loc=$windowHi  spac=$windowDiv 
y.mesh  loc=$windowLo  spac=$windowDiv 

#  Emitter 

y.mesh  loc=$emitterLo  spac=$emitterDiv 

#  Base 
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y.mesh  loc=$baseMid  spac=$baseDiv 

#  BSF 

y.mesh  loc=$bsfHi  spac=$bsfDiv 
y.mesh  loc=$bsfLo  spac=$bsfDiv 

##  Regions 

#  Cap 

#region  num=8  mater ial=Vacuum  x ,min=-$capWidth  x . max=$capWidth  y.min=$contHi 
y ,max=$contLo 

region  num=l  material=InGaN  x . min=-$capWidth  x . max=$capWidth  y.min=$capHi 
y.max=$capLo  x.comp=0.78 

region  num=2  material=Vacuum  x ,min=-$cellWidthHalf  x ,max=-$capWidth 
y ,min=$contHi  y.max=$capLo 

region  num=3  material=Vacuum  x ,min=$capWidth  x.max=$cellWidthHalf  y.min=$contHi 
y ,max=$capLo 

#  Window  [for  Ge  cell,  use  AlGaAs  with  x.comp=0.2] 

region  num=4  mater ial=InGaN  x ,min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowLo  x.comp=0.78 

#region  num=4  material=AlGaAs  x.comp=0.2  x.min=-$cellWidthHalf 

x. max=$cellWidthHalf  y.min=$windowHi  y.max=$windowLo 

#  Emitter 

region  num=5  mater ial=InGaN  x ,min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y . min=$emitterHi  y.max=$emitterLo  x.comp=0.78 

#  Base 

region  num=6  mater ial=InGaN  x ,min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$baseHi  y ,max=$baseLo  x.comp=0.78 

#  BSF 

region  num=7  mater ial=InGaN  x ,min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y.min=$bsfHi  y.max=$bsfLo  x.comp=0.78 

##  Electrodes  [for  InGaP  cell,  add  cathode  (gold)  and  remove  cathode 

(conductor) ] 

electrode  name=cathode  material=Gold  x  ,min=-$capWidth  x . max=$capWidth 
y ,min=$contHi  y ,max=$contLo 

#electrode  name=cathode  x ,min=-$cellWidthHalf  x ,max=$cellWidthHalf 

y ,min=$windowHi  y ,max=$windowHi 

electrode  name=anode  x.min=-$cellWidthHalf  x ,max=$cellWidthHalf  y.min=$bsfLo 
y ,max=$bsfLo 

##  Doping  [for  InGaP  cell,  uncomment  cap  doping] 

#  Cap 

doping  uniform  region=l  n.type  conc=$capDop 

#  Window 

doping  uniform  region=4  n.type  conc=$winDop 

#  Emitter 

doping  uniform  region=5  n.type  conc=$emitDop 

#  Base 

doping  uniform  region=6  p.type  conc=$baseDop 

#  BSF 

doping  uniform  region=7  p.type  conc=$bsfDop 
##  Material  properties 

#  Opaque  contact  [comment  out  for  InGaP  cell] 

♦material  region=8  real . index=l . 2  imag . index=l . 8 

#  Vacuum  (for  zero  reflection) 

material  region=2  index . f ile=VacuumInGa22N . opt 
material  region=3  index . f ile=VacuumInGa22N . opt 

#InGaN 

material  material=InGaN  EG300=1.1076  index . file=InGa22N . opt 
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#  Gold 

material  material=Gold  real . index=l . 2  imag . index=l . 8 


##  Models  [InGaP  cell,  1;  GaAs  cell,  5&6;  InGaNAs  cell,  7] 

models  region=l  CONMOB 

MODELS  CONMOB  FLDMOB  SRH  OPTR  PRINT 

##  Light  beams  [GaAs  bl , 0 . 55-0 . 75 , 2 00  b2 , 0 . 75-0 . 8 8 , 65 ]  0.12-2.7,50  [630,825] 
beam  num=l  x.origin=0  y . origin=$lightY  angle=90  back.refl 

power . f ile=Post Junction3 . spec . txt  \ 

wavel . start=0 . 12  wavel . end=2 . 4  wavel.num=50 

struct  outf ile=SingleCell_webf . str 
#tonyplot  SingleCell_webf . str 

solve  init 

method  gummel  newton  maxtraps=10  itlimit=25 
solve  bl=0 . 9 

##  Getting  Isc  for  I-V  curve  points 
method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 


extract  name="isc"  max 

( i . "cathode 

set 

isc=$isc*$width3d 

set 

il=$isc/10 

set 

i2=$il+$isc/10 

set 

i3=$i2+$isc/10 

set 

i4=$i3+$isc/10 

set 

i5=$i4+$isc/10 

set 

i6=$i5+$isc/20 

set 

i7=$i6+$isc/20 

set 

i8=$i7+$isc/20 

set 

i9=$i8+$isc/20 

set 

H0=$i9  +  $isc/20 

set 

ill=$il0+$isc/40 

set 

il2=$ill+$isc/40 

set 

H3=$il2  +  $isc/40 

set 

H4=$il3  +  $isc/40 

set 

H5=$il4  +  $isc/40 

set 

H6=$il5  +  $isc/80 

set 

H7=$il6  +  $isc/8  0 

set 

H8=$il7  +  $isc/8  0 

set 

H9=$il8  +  $isc/80 

set 

i20=$il9+$isc/80 

set 

i21=$i20+$isc/80 

set 

i22=$i21+$isc/80 

set 

i23=$i22+$isc/80 

set 

i24=$i23+$isc/80 

set 

i25=$i24+$isc/80-0 

.  00001 

## 

log 

outfile=In78Ga22N. 

log 

method  newton  maxtraps=10  itlimit=100 
solve  bl=0 . 95 

contact  name=anode  current 

method  newton  maxtraps=10  itlimit=100 

##  Pmax  points  [InGaP  18-25;  GaAs  15-25;  InGaNAs  13-25;  Ge  11-25] 
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solve 

ianode= 

-$i25 

bl=0 

.  95 

solve 

ianode= 

-$i2  4 

bl=0 

.  95 

solve 

ianode= 

-$i23 

bl=0 

.  95 

solve 

ianode= 

-$i22 

bl=0 

.  95 

solve 

ianode= 

-$i2 1 

bl=0 

.  95 

solve 

ianode= 

-$i20 

bl=0 

.  95 

solve 

ianode= 

-$il  9 

bl=0 

.  95 

solve 

ianode= 

-$il8 

bl=0 

.  95 

solve 

ianode= 

-$il7 

bl=0 

.  95 

solve 

ianode= 

-$il  6 

bl=0 

.  95 

solve 

ianode= 

-$il5 

bl=0 

.  95 

solve 

ianode= 

-$il4 

bl=0 

.  95 

solve 

ianode= 

-$il3 

bl=0 

.  95 

solve 

ianode= 

-$il2 

bl=0 

.  95 

solve 

ianode= 

-$ill 

bl=0 

.  95 

solve 

ianode= 

-$ilO 

bl=0 

.  95 

solve 

ianode= 

-$19 

bl=0 . 

95 

solve 

ianode= 

-$i8 

bl  =  0  . 

95 

solve 

ianode= 

-$17 

bl  =  0 . 

95 

solve 

ianode= 

-$i6 

bl  =  0 . 

95 

solve 

ianode= 

-$15 

bl  =  0 . 

95 

solve 

ianode= 

-$14 

bl  =  0 . 

95 

solve 

ianode= 

-$13 

bl  =  0 . 

95 

solve 

ianode= 

-$12 

bl  =  0 . 

95 

solve 

ianode= 

-$il 

bl  =  0 . 

95 

## 

solve  ianode=0  bl=0.95 
log  off 

extract  name="iv"  curve  (v . "anode" ,  i. "cathode")  outf ile=" IVcurveIn7 8Ga22N . dat " 
tonyplot  IVcurveIn78Ga22N.dat 

log  outf ile=doneIn7 8Ga22N . log 
log  off 
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APPENDIX  B:  MATLAB  CODE 


The  Matlab  code  provided  in  this  appendix  provided 
auxiliary  support  in  the  interpretation  and  display  Sivaco 
Atlas  log  files. 

Original  code  for  some  of  the  Matlab  functions  and 
scripts  was  obtained  from  Michalopolous  [1],  Bates,  [2], 
Green  [3],  and  Canfield  [4] .  Modifications  to  the  code  were 
made  to  support  this  thesis. 


A.  INDIUM  GALLIUM  NITRIDE  BAND  GAP  CALCULATIONS 

%InGaN  Bandgap  calculations 

%Baldomero  Garcia 

clc; 

%Constants 

EgGaN=3 .42; 

EgInN=0 .77; 
b=l .43; 

x=[0:0. 0001:1]  ; 

%  x= . 67 ; 

%  InConcentration=x 
%  GaConcentration=l-x 

%Formula  for  finding  EglnGaN 

EgInGaN= ( (1-x) . *EgGaN) + (x. *EgInN) - (b . *x . * (1-x) ) 

%Formula  for  finding  x 

solve ( ' ( ( ( 1-x) *EgGaN) + (x*EgInN) - (b*x* (1-x) ) ) -EglnGaN ' , ' x ' ) ; 

%  The  other  EglnGaN  that  we're  solving  for  is  EgInGaN=l.l 
%  EgInGaN=2 . 095 

%  Gacomp=l- (1/2/b* (EgGaN-EglnN+b- (EgGaN A2 -2  *EgInN*EgGaN- 
2*b*EgGaN+EgInNA2-2*EgInN*b+bA2+4*b*EgInGaN) A (1/2) ) ) 

%  Incomp= (1/2/b* (EgGaN-EglnN+b- (EgGaNA2-2*EgInN*EgGaN- 
2*b*EgGaN+EgInNA2-2*EgInN*b+bA2+4*b*EgInGaN) A (1/2) ) ) 

%P1  o  t 

plot (x, EglnGaN, ' b ' ) ; 

xlabel (' Indium  concentration  (unitless  ratio) '); 
ylabel ( ' InGaN  Band  Gap  (eV) ' ) ; 

title (' Indium  concentration  vs  InGaN  Band  Gap ' ) ; 
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B.  CONVERSION  FROM  DIELECTRIC  CONSTANTS  (EPSILONS)  TO 
REFRACTION  COEFFICIENTS  (N,  K) 

Optical  data  needed  to  be  provided  n-k  format.  Some  of 
the  data  obtained  was  available  in  epsilon  form.  Therefore, 
a  conversion  needed  to  be  made.  Michalopoulos  generated  the 
function  below  [1] . 

%  E2NK  Convert  the  el,  e2  pairs  into  n,  k  pairs. 

%  [n  k]  =  E2NK (el ,  e2) 

%  (c)2001  by  P.  Michalopoulos 

function  [n,k]  =  e2nk(el,  e2) 

ap  =  (el  +  sqrt(el.A2  +  e2.A2))  /  2; 
an  =  (el  -  sqrt(el.A2  +  e2.A2))  /  2; 

app  =  ap  >=  0 ; 
anp  =  an  >=  0  ; 

err  =  (app  <  0)  &  (anp  <  0)  ; 

err  =  sum (err) ; 
if  err  ~=  0 

disp ( 'ERROR! ' ) 

end 

a  =  ap  .*  app  +  an  .*  anp; 
n  =  sqrt (a) ; 
k  =  e2  .  /  (2  *  n)  ; 

C.  CONVERSION  FROM  PHOTON  ENERGY  (EV)  TO  WAVELENGTH  (UM) 

Data  was  found  in  photon  energy  and  wavelength  form. 
In  order  to  make  comparisons,  conversion  from  one  form  to 
another  was  necessary.  Functions  by  Michalopoulos  [1]  are 
given  below. 

%  EV2UM  Converts  photon  energy  (eV)  into  wavelength  (um) . 

%  (c)2001  by  P.  Michalopoulos 
function  um  =  ev2um(ev) 
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h  =  6 . 62  607  55e-34 ; 
eV  =  1 . 60218e-19; 
c  =  2 . 997  92458e8 ; 


ev  =  ev  *  eV; 
f  =  ev  /  h; 
wave 1  =  c  .  /  f ; 
um  =  wavel  /  le-6; 


%  UM2EV  Converts  photon  wavelength  (um)  into  energy  (eV) . 

%  (c)2001  by  P.  Michalopoulos 

function  ev  =  um2ev(um) 

h  =  6 . 62  607  55e-34 ; 
eV  =  1.60218e-19; 
c  =  2 . 997 9245 8e8; 

wavel  =  um  *  le-6; 
f  =  c  ./  wavel; 
ev  =  h  *  f  ; 
ev  =  ev  .  /  eV; 


D .  IV  CURVE  PLOTS  FOR  INDIUM  GALLIUM  NITRIDE  QUAD 
JUNCTION  SOLAR  CELL 


%  Code  modified  from  Canfield's  source  code 

%  Data  obtained  from  log  files  from  InGaN  quad- j unction  solar  cell 

%  Baldomero  Garcia 

clc; 

%  In20Ga80N  with  AMO  spectrum 

currentl=[0. 0176131398  0.0176031407  0.0175505986  0.0174980573 


0.0173929742 
0.0165123158 
0.0156316602 
0.0140905121 
0.0114485415 
0.0070452558 
0.0000000001]  ; 
voltagel= [0.0000000000 


0.0171728094 

0.0162921532 

0.0154114983 

0.0136501825 

0.0105678837 

0.0052839419 


0.0169526469 

0.0160719881 

0.0149711686 

0.0132098541 

0.0096872266 

0.0035226276 


0.0167324843 

0.0158518253 

0.0145308402 

0.0123291993 

0.0088065691 

0.0017613142 


0.2970913579 


1 .1822462140 


1.5075176680 


7472436590 

9875670060 

0436835120 

1013319790 

1566166780 

2006132650 


8900757790 

0047710060 

0540360290 

1132957310 

1686629450 

2108811030 


9378145180 

0194136280 

0722008890 

1239292560 

1786938920 

2190289820 


9663764170 

0322369610 

0877645130 

1419700890 

1871519430 

2257052930 


2.2313181110] ; 


In57Ga43N  with  post j unctionl  spectrum 
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current2= [ 0  . 
0.0164711457 
0.0156155026 
0.0145459485 
0.0128346588 
0.0094120828 
0.0034225761 
voltage2= [ 0 . 
0.9331338594 
0.9923675261 
1.0309568480 
1.0680958050 
1.1077533400 
1.1398643710 


0171128720  0.0171028 

0.0162572338  0 

0.0154015923  0 

0.0141181240  0 

0.0119790166  0 

0.0085564398  0 

0.0017112879 
0000000000  0.1305081 


0.9534997163 


0017892410 

0422028240 

0809744750 

1141715070 

1456000250 


783  0.01 
.0160433258 
.0151876797 
.0136903032 
.0111233713 
.0068451513 
0.000000000 
858  0.81 
. 9689573887 
.0101786220 
.0519502020 
.0915085520 
.  1247003520 
.  1505207070 


68989671  0.0166850573 

0.0158294126 
0.0149737698 
0.0132624807 
0.0102677280 
0.0051338641 

1]  ; 

92581792  0.9015131674 

0.9815965012 
1.0177449700 
1.0605082500 
1 .1002972060 
1 . 1330384420 

]  ; 


%  In68Ga32N  with  post j unction2  spectrum 


cur rent 3= [0.0141231105 


0.0137700408 

0.0130638854 

0.0123577303 

0.0109454185 

0.0084738721 

0.0042369367 


0.6842499725 

0.7483143199 

0.7789805527 

0.8157932979 

0.8511524069 

0.8818299562 


0.014 
0.0135935022 
0.0128873465 
0.0120046521 
0.0105923395 
0.0077677158 
0.0028246240 
0.131 
0.7083967768 
0.7574074704 
0.7902824668 
0.8224203640 
0.8579918718 
0.8884345077 


vo It age 3= [0.0000000000 


1131184  0.014 

0.0134169620 
0.0127108077 
0.0116515739 
0.0098861822 
0.0070615601 
0.0014123119 
1096030  0.544 

0.7248954021 
0.7653951663 
0.7999460470 
0.8337763994 
0.8639469242 
0.8940340114 


0298496  0.0139465807 

0.0132404230 
0.0125342689 
0.0112984952 
0.0091800273 
0.0056492474 
0.0000000001]  ; 

1673681  0.6320277025 

0.7377135412 
0.7725296065 
0.8083630634 
0.8431914560 
0.8738573284 
0.8988707325]  ; 


%  In78Ga22N  with  potjunction3  spectrum 


current 4= [0.0128841840 


0.0124010228 

0.0117568129 

0.0109515531 

0.0096631342 

0.0070862990 

0.0025768364 


0.5309258388 

0.5745324881 

0.6040640671 

0.6331288820 

0.6659109540 

0.6948948317 


0.012 
0.0122399704 
0.0115957622 
0.0106294488 
0.0090189240 
0.0064420895 
0.0012884183 
0.107 
0.5454971566 
0.5816962863 
0.6127739225 
0.6434973415 
0.6715037308 
0.7003322584 


vo  It age 4= [0.0000000000 


8741794  0.01272 

0.0120789198 
0.0114347108 
0.0103073439 
0.0083747191 
0.0051536706 
-0.0000000001] ; 

8032006  0.46884 

0.5569116230 
0.5881009992 
0.6203776540 
0.6521554853 
0.6808929507 
0.7050534598] ; 


75  0.0125620760 

0119178658 

0112736577 

0099852400 

0077305090 

0038652545 

43  0.5100668879 

5663876299 

5938970077 

6271079482 

6595301586 

6885228391 


312 

0. 

0. 

0. 

0. 

0. 

746 

0. 

0. 

0. 

0. 

0. 


%  Combined 

current5= [ 0  0.001991  0.004011  0.006003  0.007994  0.01001 

0.01178  0.01198  0.0126  0.01275  0.01286  0.0128841840]; 

volt age 5= [4.9793  4.9542  4.9195  4.8744  4.8185  4.7271  4.5995 

4.5723  4.4736  4.3906  4.0906  0] ; 


Pl=currentl . *voltagel  ; 
[a  b] =max ( PI ) ; 
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Pmaxl=voltagel (b) *currentl (b) ; 
Vocl=max (voltagel )  ; 

Iscl=max (currentl )  ; 

FFl=Pmaxl / (Vocl*Iscl) 

Eff 1=100* Pmaxl/ ( .1353) 

P2=current2 . *voltage2  ; 

[c  d] =max ( P2 ) ; 

Pmax2=voltage2 (d) *current2 (d)  ; 
Voc2=max (voltage2 ) ; 

Isc2=max (current2 )  ; 

FF2=Pmax2/ (Voc2*Isc2) 
Eff2=100*Pmax2/ ( . 1353) 

P3=current3 . *voltage3; 

[e  f ] =max ( P3 ) ; 

Pmax3=voltage3 (f) *current3  (f) ; 
Voc3=max (voltage3) ; 

Isc3=max (current3) ; 

FF3=Pmax3/ (Voc3*Isc3) 

Eff 3=100* Pmax3/ ( .1353) 

P4=current4 . *voltage4  ; 

[g  h] =max ( P4 ) ; 

Pmax4=voltage4 (h) *current4 (h)  ; 
Voc4=max (voltage4 )  ; 

Isc4=max (current4 ) ; 

FF4=Pmax4/ (Voc4*Isc4) 
Eff4=100*Pmax4/ ( . 1353) 

P5=current5 . *voltage5; 

[k  1] =max (P5) ; 

Pmax5=voltage5 (1) *current5  (1) ; 
Voc5=max (voltage5) ; 

Isc5=max (current5) ; 

FF5=Pmax5/ (Voc5*Isc5) 

Eff 5=100* Pmax5/ ( .1353) 

plot (voltagel , currentl , ' r ' ) ; 

%grid  on; 
hold  on; 

%plot (voltagel (b) , currentl (b) , ' o ' ) ; 
hold  on; 

plot (voltage2 , current2 , ' r ' ) ; 
hold  on; 

%plot (voltage2 (d) , current2 (d) , ' o' ) ; 
hold  on; 

plot (volt age 3 , cur rent 3 , ' r ' ) ; 
hold  on; 

%plot (volt age 3 ( f ) , current 3 ( f ) , ' o ' ) ; 
hold  on; 

plot (voltage4 , current4 , ' r ' ) ; 
hold  on; 

%plot (voltage4 (h) , current4 (h) , ' o' ) ; 
hold  on; 
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plot (volt age 5 , cur rent 5 , ' b ' ) ; 
hold  on; 

plot (volt age 5 ( 1 ) , current 5 ( 1 ) , ' or ' ) ; 
hold  off; 

xlabel ( 'Voltage  (V) ' ) 
ylabel (' Current  Density  (A/cmA2) ') 
title (' Solar  Cell  IV  Curve') 
axis ( [0  5.5  0  0.020] ) 


E.  AIR  MASS  ZERO  PLOTS 


Air  Mass  Zero  (AMO)  data  was  obtained  from  NREL  [19] 
and  plotted  using  the  Matlab  function  provided  by 
Michalopoulos  [1]  .  It  was  modified  to  provide  an  additional 
plot . 


%  AMSPECTRUMS  Plot  the  data  contained  in  spectrums .mat 
%  Check  the  source  code  and  uncomment  the  required  lines. 

%  (c)2001  by  P.  Michalopoulos 
h  =  6 . 62  607  55e-34 ; 
eV  =  1.60218e-19; 
c  =  2 . 997 9245 8e8; 
load ( ' Data\spectrums ' ) 

plot (AMO  wavel,  AMO  inti,  'b',  AM15a  wavel,  AM15a  inti,  'r'),  grid  on 
title ('AM  0  and  AM  1.5  spectra'),  xlabel (' Wavelength  [um] ' ) , 
ylabel (' Irradiance  [W  /  cmA2*um] ' ) 
axis ( [0  3  0  0.23]) 


The  code  above  was  modified  to  display  AMO  with  respect 
to  photon  energy  and  wavelength,  respectively. 

Figure  (1) 

plot (um, iradiance, ' r ' ) ; 
xlabel (' Wavelength  (um) ') 
ylabel (' Irradiance  W/cmA2*um') 
title ( ' AMO ' ) 
axis ([0  2.5  0  0.25]) 

Figure  (2) 

plot (ev, iradiance, ' b ' ) ; 
xlabel ( 'Energy  (eV) ' ) 
ylabel (' Irradiance  W/cmA2um') 
title ( ' AMO ' ) 
axis ([0  6.0  0  0.25]) 
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